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l.o  INTRODUCTION 

This  report  presents  the  results  of  the  "Rapid  Acquisition  of 
Spread  Spectrum  Study"  which  was  performed  between  October  1979  and  October 
1980  under  Army  Research  Office  contract  number  DAAG29-80-C-0005. 

1.1  Statement  of  the  Problem 

Tactical  Communications  provide  the  backbone  for  the  Command  and 
Control  functions  in  today's  military  systems;  however,  this  communications 
is  being  threatened  by  the  increased  capabilities  of  enemy  electronic  counter¬ 
measures  which  could  deny  the  availability  of  critical  communications  during  a 
stressed  situation.  The  continued  development  and  deployment  of  anti-jam  (AJ) 
systems  employing  Spread  Spectrum  communications  is  necessary  to  counter  this 
electronic  threat.  Most  spread  spectrum  systems  derive  their  spreading  signal, 
either  directly  or  indirectly,  from  a  pseudonoise  (PN)  sequence  generator.  In 
order  to  despread  the  received  signal,  and  subsequently  recover  the  transmitted 
data,  the  receiver's  PN  generator  must  be  synchronized  with  the  transmitter's 
PN  generator.  The  process  of  achieving  this  synchronization  is  commonly  called 
acquisition.  Rapid  acquisition  of  timing  synchronization  is  essential  for 
tactical  communication  systems  involving  dynamic  platforms  in  highly  flexible 
operational  scenarios.  But  addressing  spread  spectrum  acquisition  in  the 
military  environment  also  imposes  che  conflicting  constraint  that  protection 
against  intelligent  jamming  not  be  unduely  compromised  during  the  acquisition 
process. 

The  performance  of  any  acquisition  or  reacquisition  strategy  is 
a  function  of  the  timing  uncertainty,  hardware  complexity  and  code  randomness. 
Many  systems  choose  to  periodically  reduce  code  randomness  via  short  codes, 
fixed  preambles,  dedicated  acquisition  intervals,  etc.  in  order  to  achieve  a 
desired  acquisition  time.  However,  this  compromise  increases  the  vulnerability 
to  smart  jammers  which  may  make  this  approach  acceptable  only  for  acquisition 
in  a  benign  environment.  Other  means  of  reducing  acquisition  time  are  generally 
at  the  expense  of  more  costly  and/or  more  complex  hardware.  For  example,  a 
more  accurate  time  standard,  or  a  second  receiver  (such  as  GPS)  for  external 
time  synchronization,  can  be  used  to  reduce  time  uncertainty.  The  large  number 
of  variables,  some  of  which  are  scenario  based  and  some  technology  based,  makes 
it  convenient  from  both  a  concept  viewpoint  and  a  design  viewpoint  to  identify 
three  categories  of  acquisition: 
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Scenario  Dependent 


Initial  Acquisition 

•  Initial  System  Operations 

•  Initial  Network  Entry 

"Coarse"  Reacquisition 

•  Network  Entry 

t  Silent  Running 

"Fine"  Reacquisition 

•  TDMA 

•  Noise  Unlocks 

•  LOS  Outage 

Initial  acquisition  is  acquisition  in  the  truest  sense  -  initial 
establishment  of  timing  synchronization  from  total  timing  uncertainty.  In 
general  it  is  the  most  difficult  and  most  time  consuming,  and  the  one  most 
sensitive  to  scenario.  Initial  system  operations  require  acquisition  of  sys¬ 
tem  timing  as  only  one  part  of  the  overall  operation  of  establishing  communi¬ 
cations  and  must  be  accompanied  by  terminal  deployment,  antenna  pointing,  and 
other  setup  functions.  Thus,  the  initial  timing  acquisition  may  be  done  in  a 
relatively  leisure  fashion  if  it  is  not  the  pacing  item.  Also,  in  some  scenarios, 
initial  acquisition  is  performed  in  a  benign  environment  such  that  AJ  performance 
is  not  required  during  this  mission  phase  and,  therefore,  may  be  accommodated  by 
use  of  a  short  repetitive  code  or  even  a  "pre-launch"  hard  wire  connection. 

Initial  entry  of  a  platform  into  a  communication  network  may  require  a  complete 
initial  acquisition  operation  or  only  a  "coarse"  reacquisition  or  "adjustment 
of  timing"  operation  depending  upon  whether  or  not  a  priori  knowledge  of  the 
network's  timing  is  available  to  the  platform. 

Reacquisition  from  "coarse"  timing  uncertainty  is  required  at  the 
completion  of  operational  modes  referred  to  as  "silent 
running"  or  "dead  reckoning".  In  these  modes  a  vehicle  is  navigated  to  a 
target  area  by  an  onboard  system  with  the  communication  link  "silent"  and  upon 
arriving  at  the  area  the  link  must  be  rapidly  established  in  a  heavy  ECM  environ¬ 
ment.  Reacquisition  strategy  and  performance  for  these  modes  are  obviously  both 
scenario  dependent  (i.e.,  availability  of  position  estimates,  etc.)  and  tech¬ 
nology  dependent  (i.e.,  clock  stability,  reacquisition  hardware,  etc.). 

The  third  category  of  timing  acquisition  is  that  of  reacquisition 
from  "fine"  timing  uncertainty,  and  is  the  acquisition  operation  most  often 


Technology 

Dependent 


Scenario  and 
Technology  Dependent 
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performed  in  tactical  communication  systems.  It  is  performed  after  relatively 
short  periods  of  link  outage  which  occur  between  bursts  in  TDMA  systems  or  as  a 
result  of  noise  or  LOS  produced  losses-of-lock.  The  timing  uncertainty  is 
small;  rapid  reacquisition  is  required;  and  reacquisition  performance  is 

basically  only  a  function  of  receiver  technology. 


1.2  Objectives  of  the  Study 

The  technical  proposal  upon  which  the  study  was  based  identified 
four  specific  objectives;  all  four  objectives  were  achieved  during  the  study. 
The  first  objective  was  to  establish  the  performance  of  an  existing  acquisition 
technique  which  would  serve  as  a  standard  by  which  to  judge  the  techniques 
developed  during  the  study.  The  sequential  detection  acquisition  technique 
was  identified  as  the  standard  so  that  the  first  objective  was: 

Objective  1:  Establish  the  acquisition  performance  of  the 
sequential  detection  algorithm  for  a  realistic 
set  of  operating  conditions. 

The  theory  and  performance  of  the  sequential  detection  algorithm  are  presented 
in  Section  3,  along  with  those  of  the  other  serial  search  techniques.  Sequen¬ 
tial  detection  is  shown  to  provide  an  acquisition  time  advantage  of  two  to  four 
times  over  the  others  at  the  low  signal -to-noise  ratios  commonly  encountered 
in  tactical  AJ  data  links.  The  second,  and  primary,  objective  of  the  study  was 
to  identify  acquisition  techniques  which  can  provide  significantly  more  rapid 
acquisition  than  the  sequential  detection  technique,  with  a  stated  goal  of  1000 
times  improvement.  As  stated  in  the  proposal: 

Objective  2:  Postulate  and  evaluate  new  techniques  for  rapid 
acquisition  -  goal  of  X1000  improvement. 

The  matched  filter  (MF)  devices  described  in  Section  5  provide  one  key  to 
orders-of-magnitude  improvement  in  acquisition  time.  The  performances  of 
direct  search  and  MF  based  acquisition  systems  are  compared  in  Section  6  and 
an  example  presented  there  indicates  a  2000  times  advantage  for  the  MF  system. 
The  use  of  special  acquisition  protocols  or  strategies,  which  are  presented  in 
Section  8,  can  provide  many  orders-of-magnitude  improvement  in  acquisition 
time.  These  techniques  are  particulari ly  applicable  to  initial  acquisition 
or  to  reacquisition  with  extremely  large  time  uncertainty.  With  limited  time 
uncertainty  the  Code  Validity  Interval  (CVI)  methods  discussed  in  Section  9 
permit  data  demodulation  without  any  reacquisition  operation.  These  techniques 
are  quite  attractive  for  application  to  tactical  AJ  data  links  in  which  bursty 
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type  messages  are  transmitted  between  dynamic  platforms  in  a  highly  flexible 
operational  scenario. 

Another  means  of  achieving  more  rapid  acquisition  is  by  choice  of 
spreading  waveform  -  the  subject  of  the  third  objective. 

Objective  3:  Exploit  the  rapid  acquisition  potential  of  FH/DS 
signal s. 

The  acquisition  performances  of  FH  and  FH/DS  are  established  and  compared  to 
that  of  DS  in  Section  7.  The  FH  and  FH/DS  waveforms  can  be  utilized  to  reduce 
acquisition  time  by  several  orders-of-magni tude  and  thus  can  be  considered  to 
be  acquisition  aids.  Basically  FH  systems  gain  the  advantage  as  a  consequence 
of  the  fact  that  the  timing  resolution  required  in  these  systems  for  data 
demodulation  is  much  coarser  than  that  required  in  DS  systems.  Therefore, 
the  timing  uncertainty  remaining  at  the  completion  of  the  acquisition  process 
is  greater.  If  the  data  link  must  provide  a  highly  accurate  range  estimate, 
as  is  common  in  DS  applications,  then  timing  must  be  resolved  to  the  same 
degree  in  the  FH  system  and  the  DS  system,  and  the  FH  system  does  not  exhibit 
any  acquisition  time  advantage. 

As  stated  in  the  proposal  the  emphasis  was  to  be  on  rapid  acqui¬ 
sition  concepts  and  that  results  were  not  to  be  constrained  by  present  day 
hardware  limitations.  However,  a  secondary  objective  was: 

Objective  4:  Implications  to  hardware  development. 

The  key  hardware  item  for  improved  acquisition  performance  is  the  MF.  This 
is  a  rapidly  developing  technology  but  today's  state-of-the-art  devices 
(described  in  Section  5)  cannot  fully  satisfy  the  requi rements .  The  major 
deficiency  is  the  limited  integration  time  or  "length",  which  limits  the 
processing  gain. 


1.3  Organization  of  the  Report 

The  basic  approach  taken  in  this  study  to  identify  rapid  acquis¬ 
ition  techniques  has  been  to  postulate,  characterize  and  compare:  (1)  receiver 
structures  and  devices,  (2)  spreading  waveforms  and  (3)  special  acquisition 
protocols  or  strategies,  and  then  to  identify  those  characteristics  which  can 
potentially  be  exploited  to  provide  a  rapid  acquisition  capability.  The 
objective  has  been  to  identify  orders  of  magnitude  of  potential  improvement 
in  acquisition  performance  over  presently  employed  direct  search  acquisition 
methods.  At  the  same  time  the  need  to  limit  the  vulnerability  to  intelligent 
jammers  has  been  a  recognized  requirement. 
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All  of  the  concepts  presented  in  this  report  are  applicable,  in 
general,  to  all  three  categories  of  acquisition.  However,  it  is  recognized 
that  it  is  entirely  impractical  to  consider  establishing  initial  synchronization 
(starting  from  total  timing  uncertainty)  directly  from  the  fully  random  spread¬ 
ing  waveforms  commonly  employed  in  military  communication  systems.  Pragmatically, 
true  initial  acquisition  can  be  accomplished  in  an  operationally  acceptable 
amount  of  time  only  if  the  randomness  of  the  spreading  waveform  is  reduced  which 
implies  some  potential  reduction  in  AJ  protection  against  "very  smart"  jammers. 
Methods  for  achieving  this  trade-off  are  specifically  addressed  in  Sections  8 
and  9;  the  other  sections  of  the  report  are  primarily  applicable  to  the  "coarse" 
and/or  "fine"  acquisition  processes.  Section  2  presents  a  heuristic  derivation 
of  a  previously  obtained  expression  for  mean  acquisition  time.  The  advantage 
of  the  heuristic  development  is  that  it  permits  greater  insight  into  the  under¬ 
lying  interactions  between  system  parameters  such  as  detection  probability, 
false  alarm  probability,  dwell  time  and  signal-to-noise  ratio.  In  addition, 
the  problem  is  formulated  there  as  a  general  search  problem  in  order  to  emphasize 
the  broad  applicability  of  the  results.  In  Section  3  the  performances  of  the 
direct  serial  search  techniques  of  single  and  double  integration  and  sequential 
detection  are  evaluated  and  compared  under  common  operating  conditions.  The 
obvious  extension  to  parallel  search  operations  is  considered  in  Section  4. 

Ideally  one  would  expect  the  acquisition  time  to  decrease  linearly  with  the 
number  of  parallel  search  devices  employed;  however,  the  results  presented 
indicate  that  this  advantage  is  not  fully  realizable.  The  limitation  arises 
from  the  fact  that  if  all  parallel  devices  are  stepped  ahead  synchronously  then 
all  of  the  devices  are  tied  up  whenever  any  one  incurs  a  false  alarm. 

The  rapidly  emerging  technology  of  matched  filter  (MF)  devices  is 
expected  to  play  a  key  role  in  the  implementation  and  performance  of  future 
spread  spectrum  systems,  both  for  acquisition  and  data  demodulation.  Character¬ 
istics  of  typical  state-of-the-art  MF  devices  which  shown  promise  for  spread 
spectrum  applications  (Charge  Couple  Devices,  SAW  Convolvers,  Digital  Corre¬ 
lators  and  the  Hybrid  Array  Correlator  presently  under  development  at  Harris), 
are  briefly  described  in  Section  5.  Section  6  compares  MF-based  acquisition  tech¬ 
niques  to  the  direct  search  techniques  with  the  intent  of  both  illustrating  the 
potential  acquisition  gains  available  with  MFs  and  providing  some  insight  into 
the  underlying  conceptual  differences  between  the  direct  search  and  MF  techniques. 
Section  examines  the  influence  of  the  spreading  waveform  on  acquisition  performance, 
presenting  basic  acquisition  structures  and  performance  characteristics  for  direct 
sequence  (DS),  frequency  hop  (FH)  and  combined  FH/DS  spreading  waveforms.  FH  has  some 


times  been  proposed  as  an  acquisition  aid  for  DS  systems  and  this  issue  is 
thoroughly  examined.  Special  acquisition  protocols  or  strategies  are  presented 
and  analyzed  in  Section  8.  These  techniques  are  primarily  of  interest  for 
acquisition  with  extremely  large  timing  uncertainty  such  as  initial  acquisition 
and  the  trade-off  between  acquisition  time  and  randomness  in  the  spreading  wave¬ 
form  is  addressed  here.  Sections  6,  7  and  8  are  tied  toegether  in  that  Section  6 
considers  acquisition  performance  as  a  function  of  receiver  hardware,  Section  7 
as  a  function  of  spreading  waveform  and  Section  8  as  a  function  of  special 
strategies.  Section  9  describes  a  signaling  technique  which  does  not  normally 
require  any  reacquisition  operation  in  the  normal  sense  but,  as  would  be 
expected,  suffers  some  increase  in  vulnerabil ity  to  smart  jamming.  The  main 
body  of  the  report  concludes  with  a  summary  of  the  study  results,  Section  10, 
and  recommendations  for  future  work.  Section  11. 

Appendix  A  describes  an  Autogressive  Spectral  Estimation  Acquisition 
Technique  (ASEAT)  for  acquisition  of  FH  signals.  The  technique  has  a  very  rapid 
acquisition  capability  but  requires  that  a  linear  feedback  shift  register  be 
used  for  generating  the  FH  addresses.  This  makes  the  system  very  vulnerable  to 
a  relatively  unsophisticated  jammer  so  that  this  technique  is  inappropriate  for 
military  spread  spectrum  systems,  at  least  in  its  present  state  of  develooment. 
Component  codes,  described  in  Appendix  B,  offer  a  means  of  introducing  structure 
into  a  DS  spreading  signal  in  such  a  way  that  it  can  be  exploited  by  the  receiver 
to  reduce  acquisition  time.  Of  course,  this  structure  (reduced  randomness)  carries 
with  it  some  potential  reduction  in  smart  jammer  immunity,  but  these  codes  do 
represent  an  alternative  to  other  structured  waveform  methods  such  as  conventional 
short  codes.  Appendix  C  discusses  the  sources  of  timing  error  and  their  relative 
contributions  to  the  total  timing  uncertainty  during  periods  of  link  outage. 
Appendix  D  discusses  optical  processing  techniques  which  have  the  potential  for 
providing  better  MF  performance  than  the  more  conventional  MF  approaches.  These 
techniques  warrant  further  study  to  establish  their  applicability  to  military 
spread  spectrum  systems  as  acquisition  as  well  as  data  demodulation  devices. 
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MEAN  SEARCH  TIME 


An  expression  for  the  mean  time  to  acquire  a  pseudonoise  (PN)  spread 
spectrum  signal  was  developed  in  1 1  ]  using  a  Markov  chain  model.  In  this  section 
a  heuristic  approach  is  used,  resulting  in  the  same  expression,  in  order  to  ob¬ 
tain  greater  insight  into  the  acquisition  process  and  to  more  readily  identify  the 
underlying  trade-off  between  the  probability  of  acquisition  detection  and  prob¬ 
ability  of  false  alarm.  In  addition,  the  formulation  of  the  problem  is  general¬ 
ized  somewhat  to  emphasize  the  fact  that  the  result  is  applicable  to  a  larger 
class  of  search  problems  than  acquisition  of  PN  spread  spectrum  signals. 

In  the  general  context  of  a  search  through  some  uncertainty  region 

containing  q  cells,  one  of  which  is  the  correct  cell,  define 

q  =  Number  of  cells  in  the  uncertainty  region 

Pp  ®  Probability  of  detection  when  at  the  correct  location 

PpA  *  Probability  of  false  alarm  (false  detection) 

C  *  Average  number  of  cells  searched  for  detection 

Q  =  Average  number  of  false  alarms  during  the  search 

K  =  Number  of  time  units  of  false  alarm  penalty 

R  *  C  +  KQ  =  Average  total  number  of  time  units  in  the 
search 

td  *  Search  (dwell)  time  at  each  cell 
T  =  Average  search  time 

It  is  assumed  that  there  is  no  a  priori  knowledge  of  the  location 
of  the  correct  cell  within  the  uncertainty  region,  and,  at  least  initially,  that 
there  is  only  one  searcher.  The  procedure  is  to  repeatedly  search  through  the 
entire  uncertainty  region  until  the  correct  cell  is  detected.  On  each  search 
through  the  region  each  cell  is  tested  only  once  but  the  order  is  arbitrary,  in 
fact  the  order  may  be  rearranged  prior  to  each  new  search  through  the  region. 

The  mean  time  to  detect  the  location  of  the  correct  cell  is  deter¬ 
mined  as  follows.  The  probability  that  the  correct  cell  is  detected  on  the  ith 
search  through  the  uncertainty  region  is 

P [Detection  on  ith  Search]  =  PQ(1  -  PQ) 

which  is  the  geometric  probability  law.  For  this  probability  law  the  mean  number 
of  trials  required  for  detection  is  1/Pp  which  implies  that  the  correct  cell  is 
detected,  on  the  average,  on  the  (1/Pp)th  search  through  the  uncertainty  region. 
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Therefore,  on  the  average,  there  are  (1/P^  -  1)  unsuccessful  searches  through 
the  q-cell  region,  or 

^  ^  9  (2.0-1) 

cells  searched  unsuccessfully  prior  to  the  final  search  of  the  region.  The 
search  will  terminate  on  the  next  search  through  the  region  at  the  location  of 
the  correct  cell.  With  the  location  uniformly  distributed  throughout  the  region 
the  mean  location  is  at 


q  +  1 


(2.0-2) 


Therefore  the  average  total  number  of  cells  searched  is  the  sum  of  (2.0-1)  and 
(2'°'2K 


=  (2-P0)q+P0 

2  K 


(2.0-3) 


To  consider  the  effect  of  false  alarms  on  the  search  process,  note 
that  of  the  C  cells  searched  1/P^  are  correct  so  that  there  is  a  possibility  of 
a  false  alarm  on  (C  -  1/PD)  cells.  The  number  of  false  alarms,  npA,  which  occur 
in  n  =  (C  -  1/Pjj)  trials  follows  a  Binomial  probability  law. 


P 


npA  False  Alarms 
in  n  Trials 


» -  w 


n-n 


FA 


np^  =  0 ,  1 ,  2 ,  . 


,  n 


The  mean  of  the  Binomial  probability  law  is 

'  "  PFA 

so  that  the  mean  number  of  false  alarms  during  the  search  is 

Q  *  (C-^)P„ 

(2  -  P„)q  *  P0  -  2 
2  P„ 


FA 


(2.0-4) 
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For  each  false  alarm  a  penalty  of  K  time  units  is  incurred  so  that  the  average 
total  number  of  time  units  spent  in  the  search  process  is 


R  =  C  +  KQ 

*  (2  *  Po>1  +  P0  *  Kl<2  -  V"  *  PD  -  2I  PFA 

2PD 

which  can  be  manipulated  into 


R  *  2  +  (2  -  PQ)(q  -  1)(1  +  K  PpA) 

p 

0 

The  mean  search  time  is  therefore 


(2.0-5) 


R  r„ 


2  +  (2  -  PD)(q  -  1)(1  +  K  PFA) 


2  Pf 


(2.0-6) 


which  is  the  mean  acquisition  time  expression  obtained  in  [  lj. 

The  importance  of  the  heuristic  derivation  presented  here  is  that 
the  effects  of  missed  detections  and  false  alarms  can  be  readily  identified  as 
individual  contributions  to  the  mean  acquisition  time.  Consequently  the  trade¬ 
offs  between  Pp,  PFA  and  Tp  can  be  examined  at  a  more  fundamental  level.  For  example 
(2.0-3)  from  (2.0-4)  it  is  seen  that  the  total  average  number  of  cells,  C, 
searched  is  a  function  only  of  Pp  whereas  the  average  number  of  false  alarms,  Q, 
is  a  function  of  both  Pp  and  PFA.  For  Pp  very  small 


(many  searches  through  the  uncertainty  region)  and 


which  for  q  >>  1  is 


Under  these  conditions  Q  is  relatively  large  since  there  are  a  large  number  of 
opportunities  (cells  searched)  for  false  alarms.  On  the  other  hand  if  Pp  is 


2-3 


close  to  unity  then 


C 


q  ♦  i 
2 


(acquisition  on  the  first  search  through  the  uncertainty  region)  and 


Q  = 


so  that  the  number  of  false  alarms  is  relatively  small  since  there  are  not 
many  opportunities  for  false  alarms.  These  dependences  on  and  Pp^  are 
intuitively  obvious  but  it  is  convenient  to  have  them  quantatively  expressed 
for  system  design  and  analysis  purposes.  The  trade-offs  between  P^,  Pp^  and 
tp  will  be  examined  in  detail  in  the  next  section  in  the  context  of  PN  code 
acquisition. 

The  other  parameter  in  the  C,  0  and  T  expressions  is  the  number  of 
cells  in  the  uncertainty  region  q.  The  three  expressions  are  in  fact  linear 
functions  of  q,  a  fact  which  is  emphasized  by  rewriting  them  as 


c  ■( 

'2  -  pd\  ,  ♦ 
* 2 

(h) 

Q  a 

r<2  -  po>pfa1 

q  + 

['PD  -  2>PFa1 

ZPD  J 

i  J 

T 


(2  -  P„)(l  *  K  P„)td 

- re: - 


q  + 

PDrD  +  (PD  '  2^K  PFATD 

- 

L  2Pd 

Obviously  C  and  Q  and  the  resulting  mean  search  time  T  can  be  reduced  by 
reducing  q.  In  a  PN  code  driven  spread  spectrum  system  with  total  timing 
uncertainty,  q  is  the  lenqth  of  the  PN  code  or  some  multiple  of  it.  Reducinq 
the  length  of  the  code  to  reduce  acquisition  time  increases  the  vulnerability  to 
intelligent  jammers,  and  protection  against  jamming,  particularly  intelligent 
jawning,  is  the  main  objective  in  many  spread  spectrum  applications.  Another 
means  of  reducing  the  mean  search  time  is  by  defining  the  uncertainty  region  such 
that  some  a  priori  preference  is  given  to  the  location  of  the  correct  cell  within 
the  region.  However,  this  may  also  increase  the  vulnerability  to  intelligent  jam 
mers.  Basically  protection  against  intelligent  jamming  requires  uncertainty  in 
the  transmitted  signal  as  viewed  by  the  jammer  but  this  same  uncertainty  is  pre¬ 
sent  in  the  signal  as  viewed  by  an  unsynchronized  friendly  receiver.  From  the 
jawning  point  of  view  the  uncertainty  translates  into  protection;  from  the 
acquisition  point  of  view  it  translates  into  long  acquisition  time.  The  obvious 


conclusion  is  that  for  reduced  acquisition  time  the  emphasis  must  be'  on  reducing 
the  dwell  time  tq.  The  rapidly  developing  matched  filter  technologies  will  un¬ 
doubtedly  make  major  contributions  in  this  area. 

One  other  obvious  means  of  reducing  the  mean  search  time  is  by  em¬ 
ploying  more  than  one  searcher.  If  S  searchers  are  employed  then  in  theory  C, 

Q  and  T  are  reduced  by  a  factor  of  S.  However,  in  practice  this  linear  reduction 
is  usually  not  realized.  The  problem  is  that  parallel  search  operations  are 
usually  implemented  such  that  the  search  devices  are  synchronuously  stepped 
through  their  respective  portions  of  the  uncertainty  region.  Consequently  any 
time  one  search  device  suffers  a  false  alarm  all  pay  the  false  alarm  penalty.  As 
the  number  of  parallel  devices  becomes  large  the  false  alarm  penalty  becomes  a 
limiting  factor. 

2 . 1  Example:  Direct  Serial  Search 

Consider  the  single  dwell  time,  noncoherent  detector  PN  code  acqui¬ 
sition  model  shown  in  Figure  2.1-1,  After  mixing  the  received  signal  with  the  local 
PN  code  reference  the  resulting  signal  is  fiUered,  envelope  detected  and  sampled 
every  rQ  seconds  where 


The  sample  value  is  then  compared  to  a  decision  threshold.  If  a  "hit"  results 
the  system  conducts  a  verification  test  for  Ktp  seconds  which  results  in  either 
a  verified  acquisition  or  a  dismissed  false  alarm.  When  either  a  “no  hit” 
decision  results  from  the  threshold  decision  or  a  false  alarm  is  dismissed  the 
PN  code  reference  is  stepped  one  chip  or  less  and  a  new  sample  value  generated 
td  seconds  later.  This  search  procedure  continues  until  acquisition  is  achieved 
Typically  the  PN  code  reference  is  stepped  in  half-chip  increments  so  that  q  is 
twice  the  code  length.  When  the  incoming  PN  signal  and  local  reference  are  not 
in  alignment  the  input  to  the  envelope  detector  is  assumed  to  be  bandlimited 
white  Gaussian  noise  of  power  N;  when  the  PN  signals  are  aligned  the  input  is 
assumed  to  be  a  sinusoid  of  peak  amplitude  Aq  plus  bandlimited  white  Gaussian 
noise  of  power  N.  Under  these  assumptions  the  samples  are  Rayleigh  distributed 
in  the  noise  only  case 

x  _x2/2N 

Px|n  (x)  -  jf  e  x  >  0 
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(2.1-1) 


so  that  the  false  alarm  probability  is 


-Vt2/2N 


FA 


where  is  the  decision  threshold.  Expressing  in  terms  of  PpA 


VT  =  V-2N  In  PFA 


(2.1-2) 


For  the  signal  plus  noise  case  the  samples  are  Ricean  distributed  but  this 
distribution  is  difficult  to  work  with  analytically  so  that  the  approximation 

I21 


px|s+n 


(x) 


-(x-A0)c/2N 


2*A0N 


will  be  used,  with  the  result  that 

—  -(x-A0)Z/2N 


WTH 
0 


dx, 


dx 


x  >  0 


(2.1-3) 


This  approximation  to  Pg  is  quite  good  whenever  the  predetection  SNR  is  somewhat 
greater  than  unity.  In  the  following  analyses  less  accurate  approximations  will 
be  required  at  some  points  to  obtain  results  which  are  simple  enough  to  be  use¬ 
ful  in  system  design.  The  accuracy  of  these  results  will  be  tested  by  considering 
the  results  obtained  from  the  approximation  in  (2.1-3)  as  "exact"  results,  with  a 
change  of  variable  (2.1-3)  becoiugs 

V5*  -u? 

e  du 


V*1  /  (vW+1) 


_ 0_ 

^2N 


but  from  (2.1-2) 

T=  s  V-1n  PFA  (2.1-4) 

V2N  V  FA 

and  Ao/V2?T  is  the  SNR  at  the  input  to  the  envelope  detector  in  the  signal  plus 
noise  case, 

SNRed  =  i  (2.1-5) 

V'2N 
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x<» 


so  that 


ao 


PPA  -VsNRed 


_u2  (2.1-6) 
e  du 


This  expression  provides  the  trade-off  between  PQ  and  PFA  needed  to  examine  C, 

Q,  and  T,  with  SNRED  as  an  additional  system  parameter.  The  predetection  SNR  is 
in  fact  a  function  of  since 


SNR 


ED 


'  (Vo 


,Fi) 


IF 


SNR 


IF 


With  the  notation 


“  -  bif  snrif 


which  has  units  (sec*1),  then 

SNRED  =  Qtd  (2.1-7) 

Values  for  a  typical  PN  spread  spectrum  system  may  be 


SNRIF  =  -20  dB 
BIF  *  10  MHz 

so  that  a  typical  value  of  a  is 

“  =  105  sec-1 

Using  (2.1-7)  the  mean  acquisition  time  normalized  by  a  to  a  unitless  quantity, 
can  be  expressed  in  terms  of  the  predetection  SNR  as 


®T  =  2  +  (2  -  PD)(q  -  1)(1  +  K  PpA) 


This  normalized  mean  acquisition  time  is  plotted  in  Figure  2.1-2  as  a  function  ot 
SNRed  for  several  values  of  PpA  and  with 

q  =  1000 
K  =  5 


It  is  seen  that  at  lower  values  of  SNRpg  acquisition  time  can  be  reduced  by 
allowing  a  higher  PpA  whereas  at  higher  values  of  SNRpo  acquisition  time  is 
somewhat  insensitive  to  PpA>  In  all  cases  an  optimum  value  of  SNRpg  exists 
(with  a  corresponding  optimum  value  of  Tg),  with  performance  more  sensitive 
to  SNRpp  at  lower  values  of  PpA.  Figure  2.1-3  provides  for  a  more  detailed 
examination  of  the  effects  of  dwell  time  on  the  acquisition  process.  In  the 
lower  range  of  a  small  increase  in  tg  produces  a  large  decrease  in  the 
average  number  of  cells  searched,  C,  due  to  the  corresponding  increase  in  Pg. 
The  average  number  of  false  alarms,  Q,  also  decreases  simply  because  there  are 
less  opportunities  to  make  a  false  alarm.  With  both  C  and  Q  decreasing,  the 
average  acquisition  time  decreases  rapidly  even  though  the  dwell  time  is 
increasing.  For  values  of  xQ  in  the  vicinity  of  55  the  rate  at  which  C  and 
Q  are  decreasing  with  increasing  xQ  is  reduced  (implying  that  Pp  is  approaching 
unity)  so  that  T  flattens  out.  As  xQ  increases  beyond  this  region  C  and  Q 
remain  essentially  constant  (slowly  approaching  their  minimum  values  corres¬ 
ponding  to  Pq=1)  so  that  T  increases  linearly  with  Xg. 

These  results  indicate  an  optimum  value  of  SNRpg,  with  a  corres¬ 
ponding  optimum  value  of  Xg,  and  with  T  exhibiting  a  relatively  sharp  null 
about  the  optimum  value  for  large  SRNppS.  Expressions  for  these  optimum 
values  can  be  obtained  by  defining 


so  that 


and  finding  the  value  of  x  which  minimizes 


T  =  2  M2  -  PD)(q  -  l)(l  +  K  PpA) 


(2.1-9) 


(2.1-10) 


(2.1-11) 
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keeping  in  mind  that  P^  is  a  function  of  x.  Differentiating  T  with 
x,  setting  equal  to  zero  and  defining 


yields 


fi  *  (q  -  1 ) ( 1  +  K  PFA) 


x  dPD  -2Pn+  $  P  2 

.  "HT  D  ft  D 


=  0 


x  =  x 


opt 


The  approximation  for  PQ  used  earlier  is 


*  /V5 


“(u-A  )  /2N 
e  0  dx 


A  less  accurate  approximation  which  will  yield  simplier  results  is 

approximation  to  (2.1-14)  [2] 

'°°  1  -(u-A  )2/2N 

e  0  du 


which  can  be  further  approximated  by  [  3} 

p  *  4-  JL. 


In  terms  of  the  optimization  parameter  x  this  becomes 


1 


1 


-(/TrT 


PFA  '  x) 


PFA  'x 


Using  (2.1-17)  in  (2.1-13)  and  defining 


=  PFA  -  X 


yields 


FA  opt 


4Y3  -  4Y2  +  2Y  +  P 


(0+lJv^T 


FA 


For  specific  values  of  PpA  and 

0  -  (q  -  1 ) ( 1  +  K  PFA) 


respect  to 

(2.1-12) 

(2.1-13) 

(2.1-14) 

the  Gaussian 

(2.1-15) 

(2.1-16) 

(2.1-17) 

(2.1-18) 

(2.1-19) 
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numerical  techniques  can  be  employed  to  solve  (2.1-19)  for  Y  and  that 

/  - -  ,2 


ED  opt 


(V^ 


(2.1-20) 


D  opt 


ED  opt 


(2.1-21) 


The  optimum  SNRpD  as  determined  from  (2.1-19)  (labeled  ‘'approximation")  is  compared 
in  Figure  2.1-4  with  the  optimum  value  as  determined  by  iterative  evaluations  of 
(2.1-8)  using  the  more  exact  approximation  for  given  in  (2.1-6)  (labeled 
"actual").  It  is  seen  that  the  "approximation"  somewhat  over  estimates  the 
"actual"  value  at  the  higher  values  of  Pp^.  But  as  noted  earlier  the  minimum 
is  quite  broad  for  the  higher  values  of  PpA  so  that  the  "approximate"  optimum 
SNRpD  should  provide  a  sufficiently  accurate  estimate  for  system  design. 

The  foregoing  considered  the  minimization  of  mean  acquisition  time 
by  selection  of  the  predetection  SNR  or,  correspondingly,  dwell  time.  The  sys¬ 
tem  designer  may  also  have  the  option  of  selecting  the  PpA  so  as  to  minimize  T. 

This  requires  the  minimization  of  T  with  respect  to  SNR^g  and  PpA.  With 


-In  Pr 


(2.1-22) 


and,  as  before. 


differentiating  T  with  respect  to  X  and  setting  equal  to  zero,  and  with  respect 
to  Z  and  setting  equal  to  zero  yields 


v  _ 5.  -  2  P  +  ft 

x  ax  c  pd  e+i 


X  =  Xc 

z  =  z. 


(2.1-23) 


-q  +  (q  -  1)K 

aT 


FA  -  Pe 


1)K  Pr 


(2.1-24) 


X  =  XC 

z  =  z. 
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where 


aPD  -  1 


-  2 


-(X-Z)‘ 


TT  2-/T  l 

aPFA  «  0 

TT 

*p0s  _lLl  J  e'(x-z)5 

ar  2VT  l2'x  i 


ap 


-Z 


and  P  is  now 


FA  g  -2  Z  e 

az 


-Zc 

3  *  (q  -  1 ) (1  +  K  e  ) 


(2.1 -25a ) 

( 2 . 1 -25b ) 
(2.1-25c) 

(2. l-25d) 

(2.1-26) 


The  resulting  pair  of  equations  which  must  be  solved  simultaneously  to  obtain 

>2 


Xopt  and  Zopt  is 


P 


4(Z-X)J  -  4(Z-X)‘  +  2(Z-X)  +  ((J+1)vr-  e 


(Z-X)‘ 


+  2Z  =  0 
(2. l-27a) 


and 


(™t2) 


+r  2x  +  l 

e  +  (q  -  1 )K  Z-X 


+  _ -2z, 

2VT  (Z-X)' 


e"(Z’X)  =  o  (2.1 -27b ) 


Once  X_*  and  Z  .  have  been  found 
opt  opt 

SNRED  opt  =  Xopt2 


td  opt  *  !HhL0£t 


and 


FA  opt 


‘opt 


(2.1-28) 

(2.1-29) 

(2.1-30) 
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This  procedure  does  not  seem  to  be  sufficiently  simple  to  be  useful  as  a  design 
tool.  In  addition,  it  is  noted  from  the  previous  examples  (Figures  2.1-2  and 
2.1-3)  that  the  optimum  value  of  PpA  is  greater  than  0.1  (for  a  false  alarm 
penalty  of  K=5)  which  is  probably  unacceptable  from  other  system  considerations. 
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3.0  SERIAL  SEARCH  TECHNIQUES 

One  technique  used  to  obtain  coarse  code  synchronization  is  the 
serial  search,  fixed  dwell  time  acquisition  system.  This  model  assumes  there 
are  q  cells  to  be  searched,  where  q  is  related  to  the  code  length.  If  the  code 

search  is  done  in  h  chip  steps,  which  is  normally  the  case,  then  q  will  equal 

twice  the  code  length.  It  is  assumed  that  the  correct  code  phase  is  equally 
likely  to  be  found  at  any  given  position.  If  the  threshold  is  exceeded  at  any 
location,  a  verify  mode  is  entered  where  the  dwell  time  is  increased  to  gain 
confidence  that  the  correct  code  phase  has  been  located.  For  each  threshold 
crossing  there  is  the  probability  that  the  true  code  phase  has  not  been  found 

which  is  given  by  the  probability  of  false  alarm  (PFA).  This  increases  the 

average  acquisition  time  and  is  modeled  as  a  constant  times  the  fixed  dwell 
time  (KxD).  If  it  is  determined  that  the  correct  code  phase  has  been  found 
then  the  tracking  process  can  begin. 

One  method  of  characterizing  the  "goodness"  of  an  acquisition  tech¬ 
nique  is  the  time  it  takes  to  search  the  time  uncertainty  and  locate  the  true 
code  phase.  This  is  complicated  by  the  fact  that  this  time  is  a  random  variable 
and  the  statistics  required  to  fully  characterize  it  are  very  difficult  to 
generate.  A  partial  characterization,  however,  can  be  had,  by  determining 
the  mean  time  to  acquire  and  variance  associated  with  this  time. 

3. 1  Single  Integration 

Figure  3.1-1  is  a  block  diagram  of  the  single  integration  model  being 
considered.  The  input  PN  signal  plus  noise  is  multiplied  by  the  receiver's 
\  best  guess  as  to  what  the  transmitter's  PN  code  phase  actually  is.  The  resul¬ 
tant  output  signal  is  then  bandpass  filtered  (BPF)  and  applied  to  a  square  law 
device.  The  bandwidth  of  the  BPF  will  be  determined  by  whether  or  not  it  has 
to  pass  data  during  the  acquisition  process.  Also  a  different  filter  may  be 
switched  in  once  acquisition  has  taken  place.  The  squared  signal  is  then 
integrated  for  a  period  of  time  determined  by  the  probability  of  detection 
requirement,  PFA  and  the  IF  signal  to  noise  ratio.  This  value  is  then  compared 
to  a  threshold  to  determine  if  possibly  the  correct  code  phase  has  been  found. 

If  the  threshold  is  exceeded,  the  integration  process  is  continued  for  a  time 
determined  by  just  how  much  confidence  is  needed  to  initiate  the  tracking  pro¬ 
cess.  If  the  threshold  is  not  exceeded,  then  a  signal  is  sent  to  the  PN 
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generator  to  cause  it  to  change  its  phase  by  some  amount,  typically  H  chip. 
The  process  is  repeated  until  the  correct  code  phase  is  found. 


Figure  3.1-1. 

Noncoherent  Fixed  Dwell  Time  Acquisition  Model 


The  mean  acquisition  time  and  variance  for  this  system  is  given  by  the 
equations: W 


T 

aZ 


(2-PP)  (WFA> 

2P0 

r02  (!♦»„)  V<17  -  ♦  0 


(3.1-1) 

(3.1-2) 


These  equations  assume 

•  No  doppler 

•  Correlation  error  fixed  at  H  the  update  size  (Pn  is  time 

invariant)  u 

PD  *  probability  of  a  detection 
PpA  =  probability  of  a  false  alarm 
K  *  penality  paid  for  false  alarm 
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q  =  number  of  cells  to  search  (q  will  equal  twice  the  code 
length  for  one-half  chip  updates) 

=  dwell  time 

If  doppler  is  present  the  following  equations  can  be  used. 


T  =  (2  -  Pn)  (1  +  KPfa)  (Nrn) 


2Pr 


'AT, 


+  TpAfc  +  KTDAfcPpA| 


(3.1-3) 


° 2  -  (1  +  KP„)  V  OF  -  p^  ♦  p7)td2 

(^c  +  Vfc  +  KTDAfcPFA) 


(3.1-4) 


N 

^C 

Tc 

Afc 

AfCTr 


AfcKrD  = 


number  of  chips  to  be  searched 

step  size  of  search  in  fractions  of  a  chip  (typically  h) 
code  doppler  in  chips 

PN  code  phase  timing  shift  due  to  code  doppler  during 
dwell  time 

code  phase  shift  during  hit  verification 


PpA  =  probability  of  a  false  alarm 

PD  =  probability  of  a  detection 


Note:  The  algebraic  sign  of  the  code  doppler  must  be  assigned  so  that  when 
the  search  is  speeded  up  by  the  code  doppler,  Afc,  is  positive  and 
vice  versa. 
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For  a  fixed  dwell  perfect  integrator  system,  the  probability  of  detection 
and  false  alarm  are  given  by:^ 


PD  =  erfc 


B-SNR/B.f  td1 


/2SNR+1  J 


(3.1-5) 


B  =  erfc'1  (PpA) 


where 


erfc(x) 


exp  (-y  /2)dy 


(3.1-6) 


SNR  =  pre-detection  signal-to-noise  ratio 
BiF  =  pre-detection  bandwidth 

The  SNR  can  be  computed  from  the  equation: 

SNR  *  10  log  +  10  log  ^BT~T~J  '  10  1og  Rs  (3-1-7) 

C/Nq  =  carrier-to-noise  ratio 
Ts  =  1/R$  symbol  time 
R$  =  coded  symbol  rate 

The  following  is  an  example  of  the  single  integration  calculations. 
The  data  used  in  this  calculation  will  also  be  used  in  future  calculations  for 
different  acquisition  architectures  so  that  a  relative  comparison  can  be  made 
between  their  performances.  The  parameters  were  chosen  arbitrarily. 
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Assume:  Bif  =  200  KHz 

SNR  =  -9.5  dB 
P0  -  .26 
PFA  =  .0015 

From  the  relation  B=erfc~*  (PFA) 


Code  =  2047  chips 
Search  in  H  chip  steps 
Assume  no  doppler 
K  =  5 


erfc(x)  = 


1 

JTn 


03  /  2  \ 

f  exp(^f) 


dy 


B  =  2.96 


Now  using  the  equation 


P 


D 


=  erfc 


0-SNR/B  j^  t q 
.  /1+2SNR 


For  Pp  =  .26  the  argument  of  the  erfc  equals  .65. 


.65  = 


2. 96-. 112  v400x103~td 
/UZCTLTj 


Solving  for  tq  =  2  ms 

(2-P  )(1+KPfa) 

Acquisition  time  T  =  - q  xn 

2P0 

T  =  l2"-26J-j •  OOJS))  (4094)  (2xio'3) 

T  =  27.6s  mean  acquisition  time  for  SNR  =  -9.5  dB 

Figures  3.1-2,  3.1-3  and  3.1-4  were  plotted  assuming  Gaussian 
statistics.  The  probability  density  at  the  output  of  a  square  law  detector 
is  "XI  SQUARED."  If  the  detector  is  followed  by  a  filter  whose  bandwidth  is 
much  smaller  than  the  predetection  bandwidth  then  the  output  is  usually 
assumed  to  be  Gaussian. 
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Figure  3.1-2. 
Acquisition  Performance 


“  ‘  SNR  (dB) 

Figure  3.1-3. 
Acquisition  Performance 


(10‘°  see) 


Figure  3.1-4. 
Acquisition  Performance 


Figures  3.1-2  and  3.1-3  show  how  acquisition  time  varies  with 
signal  to  noise  ratio.  These  curves  were  plotted  keeping  probability  of 
detection  and  probability  of  false  alarm  constant.  A  new  dwell  time  was 
calculated  each  time  the  SNR  was  changed.  Notice  that  the  curves  cross 
over  around  a  SNR  of  0  dB.  This  shows  that  at  high  SNR  the  acquisition 
time  actually  decreases  with  increasing  probability  of  detection  and  is  due 
to  the  fact  that  the  probability  of  detection  on  the  first  pass  increases 
significantly.  For  low  SNR  values  the  probability  of  detection  on  the  first 
pass  is  quite  low  and  to  achieve  a  high  probability  of  detection  requires 
several  pulses  to  be  added  together  to  increase  the  SNR  to  a  value  where 
detection  can  occur. 

Figure  3.1-4  shows  how  the  dwell  time  increases  for  increasing 
probability  of  detection  and  decreasing  SNR. 

In  the  derivation  of  the  equation  for  acquisition  time,  the  total 
time  was  shown  to  consist  of  two  components;  one  due  to  the  time  required  to 
search  the  time  uncertainty  to  build  up  the  SNR  to  a  value  where  detection 
can  be  declared  and  one  due  to  false  alarms.  Figure  3.1-5  is  a  plot  of  the 
false  alarm  contribution  to  the  acquisition  time.  For  the  probabilities  of 
false  alarm  assumed,  the  contribution  is  very  small. 
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One  disadvantage  of  the  single  integration  method  is  that  the 
dwell  time  is  fixed  so  that  just  as  much  time  is  spent  evaluating  a  noise 
only  position  as  is  spent  for  a  potential  alignment  position.  The  double 
integration  method,  which  will  be  considered  next,  cuts  down  on  the  time 
spent  evaluating  a  noise  only  position.  This  is  achieved  by  setting  the 
threshold  and  integration  time  such  that  each  phase  is  looked  at  for  a 
smaller  period  of  time  and  the  number  of  false  alarms  increases.  The  noise 
only  positions  are  then  dismissed  quicker.  When  a  "hit"  is  detected,  a 
second,  longer,  integration  is  started  to  determine  if  the  code  phase  has 
been  found  or  is  it  just  a  false  alarm.  If  the  second  threshold  is  exceeded, 
then  tracking  may  be  initiated  or  a  third  integration  may  be  attempted  to 
verify  further  that  the  true  code  phase  has  been  found.  The  basic  idea  in 
the  above  strategy  is  to  apportion  some  of  the  false  alarm  protection  in  the 
first  integration  and  place  the  remaining  (usually  greater)  protection  in  the 
second  integration.  It  will  be  seen  that  this  double  integration  method  reduces 
acquisition  time  over  the  single  integration  technique. 

The  noncoherent  double  dwell  time  block  diagram  is  shown  in  Figure 
3.2-1.  The  following  are  the  equations  and  definitions  for  this  configuration. W 


Figure  3.2-1. 

Noncoherent  Double  Dwell  Time  Acquisition  Model 
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For  fixed  dwell  perfect  integrator  system,  the  probability  of  detection  and 
false  alarm  are  given  by : ^ 


=  erfc  B  -  SNRvBif  rD^ 
J2SH  R  +  1 


erfc"1 (Pc*) 


(3.2-5) 


Where 


erfc(x)  = 


f  yfe 


exp(-y2/2)dy 


(3.2-6) 


SNR  =  pre-detection  signal -to-noise  ratio 
=  pre-detection  bandwidth 


The  SNR  can  be  computed  from  the  equation: 


SNR  =  10  log  i  +  10  log  (g~r^|  -  10  log  Rs 


(3.2-7) 


The  following  is  an  example  of  the  double  integration  calculations. 


Assume:  B.^  =  200  KHz 

SNR  =  -9.5  dB 

PDj  =  .28 

P°2=  .9 
Pfai  =  -03 
PFA2  =  ,05 


Code  =  2047  chips 
Search  in  H  chip  steps 
Assume  no  doppler 
K  =  5 


From  the  relation  6j  =  erfc-1  (PpA1) 


Bj  =  1.88 


*  -*s+ 


Solving  for  tq1  =  615  ys 

Solving  for  ^  from  ^  =  erfc-1  ( PpA2 ^ 

|32  =  1.65 

Solve  for  t p2  as  above 

For  PD2  =  *9  the  argument  of  the  erfc  =  -1.29 


tD2  =  3.75  ms 
2-P 

T  =  "2^  [jDl+TD2  PFA1(1+KPFA2)]<’  =  [•  615xl0-3+3. 75xl0"3(  .03)  (l+5( .  05) )]  4094 

T  =  10.4s  mean  acquisition  time  SNR  =  -9.5 

The  same  basic  data  that  was  used  to  evaluate  single  integration  was 
used  to  determine  the  performance  of  the  double  integration  system.  Figures 
3.2-2  and  3.2-3  are  plots  of  acquisition  time  versus  SNR  for  various  probabilities 
of  detection.  They  show  that  acquisition  time  is  not  very  sensitive  to  overall 
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Acquisition  Time  (sec) 


Figure  3.2-2. 
Acquisition  Performance 


Figure  3.2-3. 
Acquisition  Performance 


probability  of  detection  when  the  second  probability  of  detection  is  very 
high.  Notice  the  increase  in  acquisition  time  for  the  case  of  and 
equals  .71.  The  total  probability  of  detection  is  .5  but  the  acquisition 
time  is  35%  higher  than  the  acquisition  time  for  the  Ppj  equal  .51  and  P  ^ 
equal  .99.  This  is  plotted  in  Figure  3.2-2. 
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Since  double  integration  decreased  the  acquisition  time,  with  the 
short-term  integrator  dismissing  the  low  noise  values  rapidly  and  the  long 
integrator  cleaning  up  the  false  alarms  on  high  noise  values,  it  would  seem 
that  maybe  three  integrators  would  improve  acquisition  time  even  more.  One 
would  be  very  short,  rapidly  dismissing  the  lowest  noise  values,  the  next 
would  be  a  little  longer  and  clean  up  the  false  alarms  due  to  moderate  noise 
values,  and  the  third  would  do  the  final  resolving  of  the  highest  noise  sam¬ 
ples.  In  fact,  why  stop  at  three  —  Why  not  four  or  five?  Better  yet,  the 
best  system  would  be  one  which  would  automatically  adjust  its  integration 
length  to  just  long  enough  to  make  a  satisfactory  decision  on  whatever  set 
of  noise  samples  it  happened  to  be  getting.  This,  in  effect,  is  what  sequen¬ 
tial  detection  does. 

The  sequential  probability  ratio  test  was  discovered  by  Abraham  Wald 
in  the  early  1940's  while  working  on  radar  detection.  If  one  is  sampling  a 
random  variable  and  attempting  to  decide  between  two  single  hypotheses,  then 
the  sequential  probability  ratio  test  is  the  optimum  test  in  the  sense  that  for 
a  given  PpA  and  PD  it  requires  the  minimum  average  number  of  samples  of  the 
random  variable  to  produce  a  decision.  In  our  case  minimizing  the  average  num¬ 
ber  of  samples  required  is  equivalent  to  minimizing  the  average  dwell  time  per 
code  phase  position,  i.e.,  the  sequential  detector  is  an  optimum  correlation 
detector  in  the  sense  that  it  produces  the  minimum  average  dwell  time  for  a 
given  level  of  performance. 

We  will  formulate  the  test  procedure  for  the  case  of  interest  and 
then  evaluate  its  performance.  Let  ri  be  the  independent  samples  of  a  random 
variable  and  consider  that  we  are  testing  the  hypothesis  of  signal  plus  noise 
present  versus  the  hypothesis  of  noise  only  present.  Under  these  conditions, 
the  sequential  test  can  be  stated  mathematically  by. 


n 

n  P(r./S+N) 
i  =  l  1 

n 

n  P(r./N) 
i=l  1 


(3.3-1) 
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If  An  >  decide  signal  present 

An  <  VL  decide  noise  only  present 

VL  .  An  <  Vu  compute  An+j  using  Equation  3.3-1  and  repeat  test 


where  P(r./S+N)  is  the  probability  density  function  of  the  random  variable 
r.j  given  signal  plus  noise,  and  Ap  is  the  familiar  likelihood  ratio. 

For  the  case  of  interest,  r.  is  a  sample  of  the  envelope  of  a  sine 
wave  plus  narrowband  Gaussian  noise.  The  density  function  for  r.  under  these 
conditions  is  the  Rice  density,  i.e.. 


P(r./S+N)  =  V-1  e 
a 


(3.3-2) 


2 

where  A  is  the  peak  amplitude  of  the  sine  wave  signal  and  o  is  the  noise  power 
in  the  narrowband  filter. 

The  density  of  r.  given  noise  is  found  by  letting  A  go  to  zero  in 
Equation  (3.3-2).  Using  Equation  (3.3-2)  in  Equation  (3.3-1/  we  formulate  the 

test  of  interest  as. 


(r,W) 


Simplifying  Equation  (3.3-3)  and  taking  the  log  of  both  sides  we  have. 


[1] 


1  nA 

n 


n 

Z 

t=l 


In  I 

o 


A2/2a2 


(3.3-4) 
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Equation  (3.3-4)  defines  the  sequential  test.  In  Equation  (3.3-4),  for  con¬ 
venience,  we  have  compared  the  log  of  X  to  the  log  of  the  thresholds.  A  block 
diagram  of  the  sequential  detector  is  shown  in  Figure  3.3-1.  It  is  perhaps  of 
interest  to  note  that  if  one  performs  this  derivation  starting  with  a  square 
law  detector  and  computes  the  likelihood  ratio  for  samples  taken  at  the  output 
of  a  square  law  detector  he  will  find  that  the  processor  shown  in  Figure  3.3-1 
is  still  optimum.  The  likelihood  ratio  processor  for  samples  taken  at  the  out 

put  of  a  square  law  detector  contains  a  square  root  device  before  the  multipli 
2 

cation  by  A/a  ,  i.e.,  it  computes  the  envelope  first. 

2  2 

In  Figure  3.3-1,  A  /2a  is  the  SNR  in  the  narrowband  filter  and 
represents  a  design  parameter,  i.e.,  it  is  the  design  value  of  SNR.  The  per¬ 
formance  of  the  sequential  detector  is  characterized  by  the  probability  of 
detection,  probability  of  false  alarm  and  average  number  of  samples  required 
for  a  decision.  Wald  has  shown  that  for  any  sequential  test 


•l!d  „  Tu 
PFA 


(3.3-5) 


(3.3-6) 


where  Pp  is  probability  of  detection  at  the  design  value  of  SNR  and  PFA  is 
probability  of  false  alarm. ^  ^ 

For  cases  of  interest  to  us  Pp  is  near  1  and  PpA  is  very  small  and 
under  these  conditions  we  have 


PFA  "  e 


(3.3-7) 


PD  =  !  -  e 


where  TL  is  a  negative  number. 


Equation  (3.3-7)  defines  the  detection  performance  of  the  circuit, 
however,  the  average  number  of  samples  required  (dwell  time)  must  also  be  known 
if  we  are  to  evaluate  the  acquisition  performance.  The  average  number  of  samples 
required  to  dismiss  with  noise  only  at  the  input  can  be  expressed  as  (again  from 
Rubin), 


M(0)  = 


u-W  ' 


n  ^1-Pp/ +  PF  In 


(pD/pF) 


E(z/noise  only) 


(3.3-8) 


where 


2  2  ri^ 

z  =  -r  /2a  +  1  nl  -4- 

0  a 


Evaluation  of  Equation  (3.3-8)  is  difficult  because  of  the  require¬ 
ment  to  compute  the  average  value  of  lnl.(x).  There  are  three  approaches  which 

0  2 

can  be  taken  for  finding  M(0).  First  one  can  expand  lnl^r^A/cr  )  in  a  power 
series  and  compute  the  average  value  of  z  by  computing  the  moments  of  the  Rice 
random  variable  r.  Figure  3.3-2  is  a  plot  of  lnl0(x)  versus  x.  From  this  figure 
it  is  apparent  that  for  the  square  law  approximation  to  be  accurate  the  SNR  must 
be  less  than  0  dB. 

The  second  approach  to  computing  the  expected  value  in  Equation  (3.3-8) 
is  by  numerical  integration.  This  is  time  consuming  on  a  computer.  In  addition, 
this  technique  gives  us  no  clue  as  to  the  Pp  at  values  of  SNR  other  than  the 
design  value. 

The  third  approach  is  to  simulate  the  sequential  detector  directly 
on  the  computer.  This  approach  also  uses  considerable  computer  time,  but  it 
allows  the  evaluation  of  a  wide  range  of  parameters  such  as  thresholds,  design 
SNR,  actual  SNR,  gain  errors,  and  external  accumulator  steps. 

This  simulation  has  been  run  for  a  wide  range  of  parameters.  Figures 
3.3-3,  3.3-4,  and  3.3-5  summarize  some  of  the  data.  Figure  3.3-3  is  a  plot  of  pro¬ 
bability  of  miss  versus  SNR  for  a  Design  SNR  of  0  dB.  The  curves  are  given  for 
various  lower  thresholds.  Figure  3.3-4  shows  the  mean-time-to-dismiss  with  no 
signal  versus  lower  threshold  for  various  Design  SNR's.  Figure  3.3-5  illustrates 
the  effect  of  the  external  accumulator  on  the  mean-time-to-dismiss  with  no  signal. 


Figure  3.3-2. 
Acquisition  Performance 


IOMR  TMRIIMOtQ.  Tt 


WM-I) 


Figure  3.3-4. 

Effect  of  Lower  Threshold  on  Mean-Time-to-Dismiss  (No  Signal) 


Choice  of  Parameters 

The  key  design  parameters  are  Design  SNR,  T^,  Ty,  and  number  of 
external  accumulator  steps.  The  Design  SNR  determines  the  parameter  A  in 
Figure  3.3-1.  This  parameter  is  chosen  to  be  equal  to  the  minimum  SNR  for 
which  acquisition  is  to  take  place.  In  this  case,  the  chosen  Design  SNR  is 
0  dB. 

After  Design  SNR  has  been  chosen,  T^  and  mean-time-to-dismiss  are 
linearly  related.  (See  Figure  3.3-4.)  T^  also  governs  probability  of  miss 
(Figure  3.3-3).  T^  is  chosen  through  a  trade-off  between  mean-time-to-dismiss 
and  probability  of  miss. 

The  parameter  T^  is  independent  of  the  other  parameters  when  T^  is 
chosen  to  give  a  low  probability  of  false  alarm  (PFA).  Increasing  T,.  decreases 


PFA,  but  it  also  increases  the  mean  time  to  declare  a  hit.  However,  this 
time  is  usually  small  when  compared  to  the  time  spent  dismissing  samples 
during  code  search. 

The  number  of  external  accumulator  steps  is  determined  by  the  speed 
of  the  sequential  detector  hardware.  This  number  is  made  as  small  as  possible 
because  it  tends  to  increase  mean-time-to-dismiss.  (Figure  3.3-5).  An  exter¬ 
nal  accumulator  is  used  only  when  the  sequential  detector  is  not  fast  enough 
to  test  the  accumulator  contents  on  every  detector  sample. 


Figure  3.3-5. 

Effect  of  External  Accumulator  on  Dismissal  Time 


The  data  just  presented  was  generated  on  a  computer  and  is  quite 
accurate.  It  is  not  feasible  to  date  to  get  this  degree  of  accuracy  by  using 
analytical  models.  Sometimes,  however,  it  is  desirable  to  get  an  estimate  of 
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system  performance  without  having  to  use  the  computer.  The  following  equations 
can  be  used  to  calculate  the  approximate  performance  of  a  sequential  detector. 
Calculations  are  included  for  the  same  system  parameters  used  in  the  simulation 
so  the  approximate  data  can  be  compared  to  the  simulated  data  to  see  how  well 
they  compare. 


Upper  threshold 

Tu  *  ili 
u  a 

B  =  probability  of  a  miss 
a  =  probability  of  a  false 

alarm 

Lower  threshold 

T,  V  B 

L  17a 

1-B  =  Pq  (probability  of 

detection) 

Operating  characteristic  function  * 

L(a)  =  Tyh  -  1 

(3.3-9) 


h  is  approximated  by  h  =  1  -  2 
? 

width  and  arf  is  the  design  SNR 

The  average  n 
is  given  by  the  equation 


7 


2 

where  a  is  the  actual  SNR  in  the  IF  band- 


The  average  number  of  samples  required  for  a  given  set  of  parameters 


ASN 


L(a)  logeB  +  [ 1 - L ( a )]  loge  ~ 


(3.3-10) 


This  equation  was  used  to  plot  Figure  3.3-6.  Comparing  this  to 
Figure  3.3-4  which  was  generated  by  the  computer  simulation  model  it  is  seen 
that  the  two  compare  quite  well  for  T^=-l  but  they  diverge  quite  rapidly  for 
lower  threshold  values. 

Figure  3.3-7  which  shows  probability  of  miss  versus  SNR  was  plotted 
using  the  equation: 

L(a)  =  .A  \  (3.3-11) 

A-Bn 
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Figure  3.3-6. 

Effect  of  Lower  Threshold  on  Mean-Time-to-Dismiss 
(No  Signal)  (Approximation) 


Figure  3.3-7. 
Acquisition  Performance 


Again,  this  is  only  an  approximation  and  can  be  compared  with  the  computer 
model  data  in  Figure  3.3-3. 

Figure  3.3-8  is  another  plot  of  ASN  versus  SNR.  Note  that  the  number 
of  samples  required  decreases  very  rapidly  for  decreasing  SNR,  i.e.,  it  dismisses 
the  no  signal  condition  quite  rapidly.  Notice  that  at  signal  levels  where  the 
ratio  of  the  actual  SNR  to  the  design  SNR  equal  .5  there  is  a  peak  in  the  ASN. 

At  values  of  SNR  around  this  point  sequential  detection  is  not  much  better  than 
the  fixed  length  test.  This  points  out  the  fact  that  the  design  SNR  should  be 
chosen  at  the  lowest  SNR  expected  in  operation. 

Figure  3.3-9  is  a  plot  of  acquisition  time  versus  SNR  for  the  same 
parameters  used  in  the  single  and  double  integration  examples.  This  will  allow 
a  comparison  of  the  acquisition  times  for  the  different  techniques.  This  equa¬ 
tion  was  plotted  by  determining  the  ASN  for  the  no  signal  condition  as  shown 
previously.  The  design  SNR  was  assumed  to  be  the  actual  value  as  SNR  was  varied. 
The  dwell  time  was  then  calculated  using  the  ASN  values.  These  dwell  times 
values  were  then  used  in  the  same  equation  used  to  calculate  the  acquisition 
time  for  single  integration.  The  equation  is  repeated  below  for  convenience. 

(2-P  )  (1=KPfa) 

T  =  - - qt  where  q  =  number  of  uncertainty 

2Pp  positions 

t  =  dwell  time 


From  the  derivation  in  Section  2.0  it  can  be  seen  that  the  terms 


(2-Pd)  ( 1+KP pA ) 
2K 


q 


relate  to  the  probability  of  detection  on  a  particular  pass  through  the  code 
and  to  the  time  penalty  paid  for  false  alarms  and  has  nothing  to  do  with  the 
detection  technique  used.  The  desirable  feature  about  sequential  detection 
is  that  the  dwell  time  can  vary  depending  on  the  SNR.  Once  this  dwell  time  is 
chosen  then  it  is  necessary  to  determine  how  many  passes  through  the  code  is 
required  and  what  the  penalties  are  for  the  false  alarms  that  occur.  The  above 
equation  takes  these  factors  into  account. 
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Figure  3.3-9. 
Acquisition  Performance 
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Sequential  Detection  Example  Calculations 


PD  -  .26 
PFA  =  .0015 


Design  SNR  =  -2  dB  =  .63 


T  =  M  =  -  i7, 

U  a  .0015  173 


In  173  =  5.2 


T  -  8  _  .  7  4  ~ .  . 

L  '  1-a  "  .9985  =  *74  1n  *74  =  “*3 


For  noise  only 


ASN  =  **  ^og  &  =  (8)(-.3)  _  2.4  _  6  samples 
-(a /T  -(.63)2  .3969 


Assume  sampling  at  B.f  =  200xl03 


T 


s 


NM 

B 


N  samples  per  position 
M  number  of  positions 
B  =  samples/sec 


,  (6)  (4,094) 
200xl03 


.122  s 


For  signal  +  noise 


— y  =  -6.5+2  =  4.50 


.35) 


^  =  173  (5.2) 


,.3 


L ( a )  -  —023) ’,1-1  ,  __  M9  .  >g8 

(173)  -(.74)*J  4. 69-. 91  3.78 


ASN  =  ( • 98 ) ( - . 3 )  +  [l-.98][5.2]  =  -.294+. 104 

-.1  .5-. 35 


-.015 


=  13  samples 


.  (13H4094)_.  _266$ 


S  200xl03 
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3.4 


Doppler  Limitations 

One  potential  problem  that  can  degrade  acquisition  performance  is 
doppler.  In  any  system  that  includes  RPVs  and  aircraft,  this  is  definitely 
a  consideration.  It  is  both  interesting  and  informative  to  evaluate  the  typical 
detection  system  investigated  thus  far  in  a  doppler  environment.  The  incoherent 
and  coherent  detectors  will  be  investigated.  Their  performance  will  then  be 
compared  to  that  of  the  convolver  in  Section  5.3.1. 


Figure  3.4-1. 
Noncoherent  Detector 


The  received  signal,  as  shown  in  Figure  3.4-1,  is 


r(t)  =  /IS  C(t)  cos  ( [aj+co^jt) 


where:  S  =  signal  power 

C(t)  =  PN  code  [±l] 
Wj  =  doppler 


(3.4-1) 
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and  the  receiver's  replica  signal  is 


S(t)  =  C(t+x)  cos  (ojt+0) 


(3.4-2) 


where:  t  =  time  offset  of  PN  replica 

t  =  carrier  phase  offset 

It  will  be  assumed  that  the  data  (if  present)  is  at  a  rate  low  enough 
so  that  it  may  be  neglected  without  consequence. 

The  question  of  concern  is  how  ujd  effects  the  permissible  value  of 
integration  time  T.  The  outputs  of  the  multipliers  are  given  by 

Xc  =  r(t)  s(t)  =  j^S  C(t)  cos  [(a)-hjJd)t]j  C(t+x)  cos  (wt+0) 


Xc  =  2S  C(t)  C(t+x)  cos  [(w+ajd)t]  cos  [(ait+0)] 
using  the  identity 

cosa  cose  =  cos(a-e)  +  %  cos(a+e)  (3.4-3) 

Xc  =  2S  C(t)  C(t+x)  j%  cos  |>dt-d)]+  J5  cos  [(2a)+u)d)t+d]| 

Xc  =  S  C(t)  C(t+x)  J^cos  [(a>dt-j5)]  +  cos  [(2w+u>d)t+0]J  (3.4-4) 


The  output  of  the  second  multiplier  is  given  by 


Xs  =  r(t)  S(t+  J) 

Xs  =  j/T?  C(t)  cos  [u+wd]tj  |  J?S  C(t+x)  sin  (tut +0) 
Xs  =  2S  C(t)  C(t+x)  cos  ([a>+(ud]t)  sin  (cut+0) 
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using  the  identity 


cosa  sin8  =  ^  sin  (a+8)  -  H  sin  (a-6) 

Xs  =  2S  C(t)  C(t+r)  j>s  sin[(2aj+wd)t+0]  -  %  sin[u)dt-d] 

Xs  =  S  C(t)  C ( t+T )  j sin  [(2u)+cad)t  +  t]-  sin[iDdt-0]|  (3.4-5) 

If  oj»1/T  we  can  invoke  the  Riemann-Lebesgue  Lerrana  to  neglect  the  integrals 
with  2o)  terms  so  that 

T/2 

Yc  =  S  J  C(t)  C ( t+T )  cos  (todt-0)dt  (3.4-6) 

-T/2 


T/2 

=  S  cos  (j5)  J  C(t)  C(t-hr)  cos  ((Ddt)dt 
-T/2 

T/2 

+  S  sin  (f5)  J  C(t)  C(t+x)  sin(wdt)dt 
-T/2 


T/2 

Y$  =  -S  J  C(t)  C ( t+T )  sin  (wdt-?S)dt  (3.4-7) 

-T/2 

T/2 

=  -S  sin  (0)  J  C(t)  C(t+x)  cos  (codt)dt 
-T/2 

T/2 

it)  J  C(t)  C ( t+T )  sin  (o)dt)dt 
-T/2 


+  S  cos 


While  it  is  not  possible  to  evaluate  the  two  integrals  in  general,  the 
results  are  certainly  just  two  constants  a  and  B  where 


a  =  J  C(t)  C(t+r)  cos  (wdt)dt 
-T/2 


(3.4-8) 


B  =  J  C(t)  C(t+t)  sin  (LJcjt)dt 
-T/2 


(3.4-9) 


For  this  assumption,  the  outputs  of  the  two  integrators  can  be  represented 
by 


Yc  =  S  a  cos  0  +  S  6  sin  0 


Ys  =  S  8  cos  0  -  S  a  sin  0 


(3.4-10) 


(3.4-11) 


The  outputs  of  the  two  channels  are  then  summed  yielding  the  following  results 


2  2 
Z  -  YCZ  +  Y/ 


(3.4-12) 


Z  =  J^S  a  cos  0+S  B  sin  0^  +  8  cos  0-S  a  sin  dj 

2  T  2  2  2  2  1 

Z  =  S  la  cos  0  +  2a3  sin  0  cos  0  +  8  sin  0J 

+  S2  cos2  0  -  2  aB  sin  0  cos  0  +  a2  sin2  0J 

Z  =  S2ja2Jcos2  0  +  sin2  0J  +  B2Jcos2  0  +  sin2  0 J  j 


Z  =  S' 


^a2  +  B2j  since  cos2  0  +  sin2  0=1 


(3.4-13) 
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If  we  now  assume  that  the  doppler,  co^,  is  zero  and  the  PN  code  period  is 
T  then 

T/2  T/2 

a  =  J  C(t)  C(t+x)  cos  w^t  dt  =  J  C(t)  C(t+r)  dt 
-T/2  -T/2 


This  is  immediately  recognized  as  the  autocorrelation  function  Rc(t). 

a  =  Rc(t)  (3.4-14) 

8  =  0 

Therefore  the  output  can  be  written  as 

Z  -  S2[o2  +  B2] 

2  2 

Z  =  S  R  c (t )  outDut  for  no  doppler  (3.4-15) 

If  it  is  now  assumed  that  there  is  doppler  but  the  codes  are  in  alignment 
then  the  output  can  be  written, 

T/2 

a  ■  J  C(t)  C ( t+T )  cos  w^t  dt 
-T/2 

for  code  alignment  t*0 
T/2 

a  =  f  C2(t)  cos  wdt  dt  but  C2(t)  =  1 
-T/2 
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T/2 

a  =  J  cos  (u^t  dt  = 
-T/2 


sin  i^t 
0)  . 


T/2 


•T/2 


cj  .T 

_  2  sin  -g-  sin  TrfdT 


a  =  T  sine  fdT 


T/2 

0  =  /  C(t)  C(t+T)  sin  Wjjt  dt 

-T/2 

for  x=0 

T/2 

3  =  J  sin  dt 
-T/2 


(3.4-16) 


since  sin  w^t  is  an  odd  function  and  the  limits  of  integration  are  symmetrical 
then  0=0.  Now  the  output  can  be  written  as  follows. 


Z  =  S2  [a2  +02] 

Z  =  S2  T2  sine2  frfT 


(3.4-17) 


This  equation  shows  how  the  output  is  affected  by  doppler.  Note  that  the  output 
goes  to  zero  for  fdT=l,  which  is  the  same  results  achieved  with  the  convolver. 
Therefore  the  convolver  has  no  advantage  over  the  conventional  quadrature  method 
of  incoherent  detection. 
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SECTION  4.0 

PARALLEL  SEARCH  OPERATION 


4-0 
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PARALLEL  SEARCH  OPERATION 


Thus  far  all  of  the  techniques  investigated  have  been  serial  search 
techniques  where  one  code  phase  at  a  time  is  checked  for  alignment.  A  natural 
extension  would  seem  to  be  to  use  two  or  more  paths  to  search  more  than  one 
code  phase  at  a  time.  Hopefully  the  acquisition  time  would  decrease  in  direct 
relation  to  the  number  of  parallel  paths  used.  This  section  investigates  the 
possibilities  of  N  parallel  searches  and  looks  at  what  the  limitations  are. 

Rayleigh  Assumption 

Assume  that  the  probability  of  dismissing  an  individual  false  align¬ 
ment  results  in  a  Rayleigh  distribution  in  time 

2 

P.  (t)  =  4-  e_t  /2X  tlO  (4.0-1) 

td  A 


The  mean  time  to  dismissal  is  t 


and  the  variance  is  a* 


=  \  (2  -  Tj-)  . 


The  cumulative  distribution  function  for  td  is  given  by 


'!  VM3/ 


-x2/2A 
e  dx 


t2/2A 

■  /  e  dE  ■ 1  - 
o 


-t2/2A 

e 


(4.0-2) 


where  S  -  x2/2A  so  that  d£  =  2xdx/2A  =  xdx/A. 

If  we  have  N  searches  progressing  in  parallel  with  a  common  control, 
then  all  N  alignments  must  be  dismissed  before  any  can  be  dismissed.  Assuming 
all  searches  have  independent  noise  the  probability  of  dismissing  all  N  searches 
by  time  t  is  the  N-fold  product  of  dismissing  each  by  time  t,  or 


i 

\ 


(4.0-3) 


The  probability  density  function  for  t^  is  given  by  the 
of  the  cumulative  distribution  function  as 


Solving  for  the  expected  value  of  tN  gives 


t 


N 


-t2/2A 

e 


-t2/2A~|N"1 

_1  -  e 


dt 


Defining 


A  ? 

5  ■  t  /2A 


so  that  d£ 


2tdt 


/2A 


gives 


t 


N 


N  71 


f 


d£ 


N-l 

=  wm  £  (-i)k 

k=0 


derivative 


(4.0-4) 
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From  Gradshteyn  and  Ryzhik,  page  317,  equation  317 


f  *h  e-|,x  dx  = 

Jo  2W  ' 


[Re  >  0] 


(4.0-6) 


Hence 


N-l 


V»#X  <-')k (V) 


k=0 


(4.0-7) 


f  tt\  _ 

Recognizing  that  y  -4-  =  t .  and  rewriting  gi 


gives 


*N  ■  'd  N  I  <-1)k'1  [Vl 


k=l 


k-1  (H-l\  1 

k3/2 


(4.0-8) 


=  tH  V  (-l)*'1  -Wtlll _  1 

(N-k):  k(k-i) :  ** 


k=l 


N 


l=h 


_td  ^  /¥ 

k=l 


k-1 


(!) 


Exponential  Assumption 

Assume  that  the  time  required  to  dismiss  a  false  position  is  an 
exponentially  distributed  random  variable 


Pt(a)  -  j- 


a  -  0 


(4.0-9) 
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with  mean 


/  t  *-“/e  * 


B  I  C  e%  d?  =  B 

*6 


r  A  a 
K  =  B 


(4.0-10) 


The  probability  that  the  position  will  be  dismissed  by  time  t  is 


p|td<to)  ■/  T  e'“/e 


[-  e-a''6]t°  =  (1  -  e'VB) 


(4.0-11) 


If  there  are  N  searches  going  on  in  parallel,  the  probability  that  they  will 


all  be  dismissed  by  time  tQ  is 


P  all  dismissed 


]■“- 


e^'V 


(4.0-12) 


Thus  the  distribution  of  tQ  is 


ptp(‘)  ■  s-  j(1  -  e't/B)N 


(4.0-13) 


=  N  (1  -  e"t/e)N_1  (-i-  e"t/e) 


-  I  e-t/e[l  -  e-W]*-1 


(4.0-13) 


Hence  the  mean  time  to  dismiss  N  positions  is 


=  (-1)N_1  NB  r(2) 


N-l 


I 

k=0 
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MATCHED  FILTER  TECHNOLOGY 


The  key  implementation  consideration  for  future  AJ  modems  will 
be  the  type  of  matched  filter  technology  employed  for  acquisition  and  data  demod¬ 
ulation.  Due  to  the  increasing  drive  towards  greater  AJ  data  link  flexibility 
(wider  range  of  data  ranges  and  vehicle  dynamics,  etc.)  versatility  of  the  matched 
filter  will  be  a  consideration  as  well  as  the  obvious  issues  of  process  gain, 
cost,  size,  etc.  The  candidate  matched  filter  technologies  are  CCD,  digital 
correlators  and  SAW  convolvers.  In  addition,  a  correlator  device  presently  under 
development  at  Harris,  consisting  of  parallel  correlators  performing  both  analog 
and  digital  functions,  offers  considerable  promise. 

The  digital  correlators,  CCDs  and  SAW  convolvers  are  generally 
considered  to  satisfy  the  low,  medium  and  high  PN  rate  requirements,  respectively. 
However,  CCDs  seem  to  be  losing  out  since  demonstrated  performance  to  date  has 
not  lived  up  to  expectations  and  the  digital  correlators  are  moving  up  into  the 
medium  PN  rate  range.  The  utilization  of  convolvers  imposes  some  additional 
hardware  requirements  in  that:  (1)  the  reference  PN  code  must  be  reversed  in  time 
since  correlation,  not  convolution,  is  the  required  operation;  (2)  the  output  of 
the  convolver  is  at  twice  the  frequency  of  the  input  and  hence  the  correlation 
pulse  is  compressed  to  half  the  normal  width.  This  compression  requires  that 
post-convolver  processing  be  at  twice  the  normal  speed,  at  least  for  a  portion 
of  the  circuitry. 

5.1  CCD 

An  example  of  the  CCD  implementation  of  a  PN  matched  filter  is 
the  Fairchild  PTF-2.  It  is  a  512-stage  CCD  delay  line  with  taps  on  511  of  the 
stages.  Basically,  the  device  generates  the  correlation  value  by  modifying  the 
sign  of  the  analog  signal  sample  at  each  tap  by  the  sign  of  the  corresponding  PN 
chip  and  summing  the  resulting  values.  This  is  accomplished  by  using  the  polarity 
of  each  PN  chip  to  switch  the  corresponding  signal  value  onto  either  a  positive 
or  negative  summing  bus  and  then  subtracting  the  two  bus  values.  The  device 
actually  contains  dual  banks  of  PN  shift  registers,  switch  banks  and  summing 
buses  as  shown  in  Figure  5.1-1.  The  CCD  delay  line  taps  are  interleaved  into  the 
two  bants  so  that  BPSK,  QPSK  or  OQPSK  signal  formats  can  be  accommodated.  The 
May  line  can  operate  at  sampling  rates  of  up  to  10  MHz  and  the  MOS  shift 
•  ;  .vo  at  up  to  2.5  MHz  each,  rates  which  are  comparable  to  those  of  the 
•  r '-relators .  The  primary  advantage  of  the  CCD  is  the  relatively  long 
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.1-1.  511-Stage  CCD  Matched  Filter 


length  available.  In  fact  the  PTF-2  Drovides  the  delayed  analog  samples  and 
PN  reference  signals  as  outputs  so  that  two  or  more  devices  may  be  cascaded  to 
form  an  even  longer  matched  filter.  However,  charge  leakage  places  a  limitation 
on  the  maximum  length  achievable. 

5.2  Digital  Correlator 

The  state-of-the-art  in  digital  correlators  is  exemplified  by 
the  TRW  TDC1023J,  a  64-bit  monolithic  correlator  whose  functional  block  diagram 
is  presented  in  Figure  5.2-1.  In  this  device  the  signal,  after  one  bit  quanti¬ 
zation,  is  serially  loaded  into  the  "A"  register  and  the  PN  reference  is  serially 
loaded  into  the  "B"  register.  Contents  of  the  "B"  register  are  then  parallel 
dumped  into  the  "R"  holding  registerand  the  correlation  between  the  "A"  and  "R" 
registers  calculated.  The  effective  length  of  the  correlator  is  adjusted  by 
loading  "0"s  into  the  “M"  register  at  locations  which  are  to  be  unused.  The 
bit-for-bit  results  of  the  active  exclusive-OR  operations  are  digitally  summed 
to  produce  a  7-bit  word  indicating  the  number  of  positions  which  are  matched. 

This  7-bit  word  is  compared  to  the  user-set  threshold  value.  Whenever  the  cor¬ 
relation  sum  equals  or  exceeds  the  threshold,  the  threshold  flag  goes  high. 

Each  of  the  three  shift  registers  may  be  indeoendently  clocked  at 
shift  speeds  up  to  40  MHz  but  the  maximum  rate  of  correlation  between  the  "A" 
and  "R"  registers  is  20  MHz. 

5.3  SAW  Convolver 

Considerable  development  efforts  have  been  devoted  in  recent  years 
to  the  utilization  of  SAW  devices  as  matched  filters  for  spread  spectrum  signals. 
Current  state-of-the-art  SAW  convolvers  exhibit  time-bandwidth  products  ranging 
from  100  to  2000  depending  on  the  material  and  fabrication  technology  employed. 
Basically,  the  convolution  operation  is  obtained  from  the  non-linear  parameters 
of  the  piezoelectric  substrate.  The  degree  of  non-linearity  and  the  efficiency 
of  the  interaction  between  the  two  acoustic  waves  can  be  enhanced  through  several 
methods.  The  air  gap  convolver  employs  a  semiconductor  device,  separated  a  small 
distance  from  the  piezoelectric  substrate,  to  enhance  the  non-linearity.  This 
approach  results  in  the  highest  performance  to  date  but  requires  a  complex  and 
costly  fabrication  process  and  exhibits  serious  mechanical  and  temperature 
instabilities.  At  this  point  in  time  it  appears  that  the  SAW  convolvers  which 
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rely  solely  on  the  non-linear  piezoelectric  elastic  parameters  of  the  sub¬ 
strate  for  the  non-linear  interaction  (the  so-called  elastic  convolvers)  are 
better  suited  to  military  applications. 

The  performance  of  the  basic  elastic  convolver  can  be  improved  by 
compressing  the  acoustic  waves  prior  to  their  reaching  the  convolution  region. 

A  simplified  diagram  of  such  a  SAW-Beam  Width  Compression  Convolver  (SAW-BWCC) 
is  shown  in  Figure  5.3-1.  Beam  width  compression  of  15:1  has  been  achieved 
with  significant  resulting  improvement  in  efficiency  and  dynamic  range.  The 
efficiency  of  the  non-linear  interaction  is  a  function  of  the  ratio  of  length 
L  to  width  W  of  the  convolution  region.  More  importantly  from  a  spread  spec¬ 
trum  processing  gain  viewpoint  is  that  the  integration  time  is  directly  pro¬ 
portional  to  the  length  L,  with  a  current  maximum  of  about  10  microseconds. 

Two  techniques  presently  used  for  beam  width  compression  are  parabolic  horn 
focusing  with  aluminum  metal  film  on  the  substrate  and  multiple  strip  couplers 
(MSC)  parallel  to  the  input  transducers.  The  parabolic  horn  approach  results 
in  somewhat  better  bandwidth  and  center  frequency  characteristics  but  also 
typically  requires  about  ten  times  greater  overall  physical  length.  Table  5.3-1 
summarizes  performance  characteristics  of  the  current  SAW  convolvers.  SAWTECH 
claims  the  capability  of  fabricating  a  SAW-BWCC  using  MSC  beam  compression  with 
300  MHz  center  frequency,  100  MHz  bandwidth  and  10  microseconds  integration  time. 

5.3.1  Effects  of  Doppler  and  Random  Phase  Variations 

The  input  to  the  convolver  is  represented  by  f(t)eJ(lllt,  while 
g(t)wJ(‘’2t  represents  the  locally  generated  reference  signal,  as  shown  in 
Figure  5. 3. 1-1.  The  waveform  at  any  time  t  and  position  x  can  be  represented 
by  [l] 


dV  =  (Ej  +  E2)(E1  +  E2)* 


(5. 3. 1-1) 


assuming  a  second  order  nonlinearity.  Ej  and  E2  are  given  by  the  real  parts 
of  the  following: 


*1  -  -  V)e 


jfwjt-kjX) 


g(t  +  J)e 


y  j  ( co2t+k2x ) 
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Beamwidth  Compression  Convolver 


Non-Linear  Medium 


Output 

Figure  5.3. 1-1. 

SAW  Convolver 

In  these  two  equations,  k^  and  k2  represent  propagation  constants  given  by: 


The  quantities  and  represent  the  input  and  reference  frequencies  to  the 
convolver  respectively.  The  velocity  of  propagation  in  the  acoustic  material 
is  given  by  v.  Performing  the  indicated  multiplication  produces  the  following 
results 


dV-y(t-J) 


"jfw.t-k.x)  -jCu.t-k.x) 
e  +e  11 


(5.3. 1-2) 


. . —  +g(t+y) 


'  j(o)2t+k?x)  -j(co?t+k?x) 
e  c  +e  c 


+  f(t>7)  9(t+J)  cos  [(w1-oJ2)t-(k1+k2)x]  +  f(t-£)  g(t+£)  cos[(u)1+W2)t-(K1-k2)x 


In  this  expression,  the  only  term  of  interest  is  the  following: 


f(t~)  g(t+J)  cos  ^(u)1+w2)t  -  (k1-k2)xj  . 

There  are  two  methods  that  can  be  used  to  detect  this  sum  frequency  term. 

If  (*>!  is  not  equal  to  u>2  then  an  interdigital  transducer  with  a  finger  pair 
spacing,  d,  is  capable  of  responding  to  a  signal  with  a  propagation  constant 

Vo>2 

d"  V 

such  that  kcjd=2ir. 

If  on  the  other  hand  equals  w2,  then  the  waveform  does  not  vary  with  x  and 
can  be  detected  between  metal  films  laid  down  on  the  top  surface  and  lower 
surface  of  the  piezoelectric  substrate.  This  is  the  configuration  of  most 
interest  in  the  acquisition  process,  where  a  reference  waveform  is  varied 
until  it  exactly  matches  the  incoming  signal. 

The  metal  film  laid  down  on  the  substrate  performs  a  spatial  inte¬ 
gration  of  the  resultant  signals.  This  yields  the  following  expression  at  the 
output  of  the  convolver. 

L 

2  x  „ 

y(t)  =  .  f(t-^)  g(t+J)  cos  u>$t  dx  (5. 3. 1-3) 

*'"7 


where 


=0),  +  0)o. 
s  1  2 


By  making  a  change  of  variable  the  detected  output  can  be  represented  by: 


Lett-V=Tf+L 

t+.2V 


y(t)  =  v 


/ 

4 


f(x)  g(2t-x)  dx 


(5.3. 1-4) 
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This  represents  the  conventional  convolutional  integral. 

Suppose  now  that  the  input  signal  is  received  having  undergone  a 
doppler  shift.  How  is  the  output  of  the  convolver  affected?  Since  ^  is  not 
equal  to  then  kd  is  not  equal  to  zero.  The  convolver  output  can  now  be 
expressed  as: 

L 

2 

y(t)  - j  f(t-^-)  g(t+^)  cosht  -  kdxjdx  .  (5. 3. 1-5) 

T 


For  purposes  of  illustration,  assume  f { t )  and  g(t)  are  10  microsecond  pulses 
that  have  undergone  a  doppler  shift.  The  above  integral  reduces  to  the 
following. 


Since  the  limits  of  integration  are  symmetrical  and  sin  k^x  is  an  odd  function, 
this  term  will  integrate  to  zero,  leaving  only  the  cos  kdx  term.  Evaluating 
the  resulting  expression,  the  convolver  output  is  given  by: 


y(t)  = 


(T.-2t)v  sin  Trfd  (T.-2t) 
wfd  (T.-2t) 


COS  0)  t 


(5.3. 1-6) 
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In  this  example  the  input  signal  is  symmetrical  so  that  the  convolu¬ 
tion  integral  also  represents  the  autocorrelation  function.  This  can  be  seen 
for  the  case  of  no  doppler  as  shown  in  Figure  5. 3. 1-1  A.  The  output  is  the 
familiar  triangle  function,  which  is  the  autocorrelation  function  for  a 
square  pulse.  Figure  5. 3. 1-2  shows  the  output  for  a  10  kHz  doppler.  The 
output  is  normalized  with  the  value  at  no  doppler.  The  basic  shape  of  the 
waveform  does  not  change  and  the  peak  amplitude  is  reduced  only  slightly. 

Figure  5. 3. 1-3  shows  the  output  for  a  50  kHz  doppler.  The  waveform  is 
becoming  rounded  which  makes  it  more  difficult  to  determine  when  the  two 
waveforms  are  in  alignment.  Note  that  there  is  also  a  significant  decrease 
in  the  amplitude  of  the  output.  As  the  doppler  is  increased  a  point  is  reached 
where  the  peak  is  no  longer  at  the  origin.  Figure  5. 3. 1-4  clearly  shows  this 
for  a  doppler  of  80  kHz.  Here  the  doppler  causes  the  convolver  output  to  peak 
at  t  equal  to  2  ys  giving  a  false  indication  of  alignment. 

The  loss  in  convolver  peak  output  versus  doppler  is  plotted  in 
Figure  5.3. 1-5.  The  3  dB  point  is  about  60  kHz  for  this  example. 

Suppose  now  that  in  addition  to  doppler  the  input  signal  contains 
a  random  phase  term.  How  is  convolver  performance  affected?  Going  back  to 
the  original  equations,  the  input  vector  E1  is  now  represented  as  shown. 


E.  =  f(t--)  e 
1  v 


j^t-cti-kjx) 


(5.3. 1-7) 


The  reference  input  g(t)  is  assumed  to  be  a  perfect  replica  of  the  input  without 
the  random  phase  term.  If  the  same  process  is  repeated  as  before  with  the  new 
Ejthen  the  sum  frequency  term  of  interest  can  be  represented  as  follows: 

L 

2 

y(t)  =  f  f(t-J)  g(t+^)  cos[wst  +  <j>  -  kdxjdx  .  (5. 3. 1-8) 

•'L 


The  cosine  term  can  be  expanded  to  simplify  the  expression. 
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Figure  b.j. 1-1A. 

Normalized  Function  for  Output  of  Convolver  with  the  Input 
and  Reference  Waveforms  Being  10  ys  Pulses 


Figure  5. 3. 1-2. 

Normalized  Correlation  Function  for  Output  of  Convolver  with  the  Input 
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and  Reference  Waveforms  Being  10  ys  Pulses  with  50  KHz  of  Doppler 


(dB) 


and  Reference  Waveforms 


y(t)  =  J  g(t+*)  cos  |^st-kdx)  cos  <f>  -sin  (wst-kdx)  sin^jdx 


Since  the  limits  of  integration  are  symmetrical  and  sin(wst-kdx)  is  an  odd 
function,  it  will  integrate  to  0.  This  leaves  the  following  expression 
for  the  output  of  the  convolver. 


cos  r  f(t-^)  g(t+J)  cosJw$t-kdx]  dx 


(5.3. 1-9) 


This  is  recognized  as  the  same  expression  that  was  developed  previously 
multiplied  by  the  cosine  of  the  random  variable.  For  the  case  of  a  square 
wave  input  with  doppler  and  random  phase,  the  output  can  be  written  immediately 


y  ( t)  = 


(T.-2t)v  sin  Trfd  (Ti-2t) 
7Tfd  (Tr2t) 


cos  4>  cos  o)st 


(5.3.1-10) 


The  expression  for  the  output  now  contains  a  random  variable.  Therefore  a 
meaningful  description  of  the  output  is  given  by  the  mean  and  mean  square  value 
The  mean  value  is  given  as  follows;  assuming  the  random  variable  <f>  is  uniformly 
distributed  between  0  and  2tt. 


(T.-2t) v  sin  irf .  (T.-2t) 

E  y(t)  =  E  -  cos  <()  cosco  „t 

Trfd  (T.-2t)  s 


(Ti--2t)v  sin  itfd  (T..-2t) 
TTf^  (T,-2t) 


cos  4>  cos  u)$t 


E[y(t)]  =  0. 


The  mean  value  of  the  output  is  zero.  Now  determine  the  mean  square 
value.  This  is  given  as  follows: 


:[y2(t)]  = 


(T.-2t)v  sin  TTfd  (T.-2t) 
rrfd  (Ti-2t) 


COS  u>st 


2tt 

\  r  2 
2^  J  cos  ?  d£ 


The  cos2(j>  term  can  be  expressed  as  . 

The  cos2<j)  term  will  integrate  to  zero  leaving  the  following  expres¬ 
sion  for  the  mean  square  value  of  the  output. 

i2 


[y2(t)]  = 


(T-j -2t ) v  sin  TTfd  (Ti - 
~~  TTfd  (Tr2t) 


i-2t)T 


2  . 

COS 


(5.3.1-12) 


This  expression  is  evaluated  for  various  values  of  doppler  as  shown  in  Figure 
5. 3. 1-6  thru  Figure  5.3. 1-9.  The  affects  caused  by  doppler  and  random  phase 
are  similar  to  those  caused  by  doppler  alone.  The  peak  output  decreases  and 
the  waveform  gradually  changes  shape  until  the  peak  is  no  longer  at  the  origin 
but  moves  to  t  equals  2  ys  as  shown  in  Figure  5. 3. 1-9.  Figure  5.3.1-10  is  a 
plot  of  the  loss  in  convolver  mean  square  output  due  to  doppler  and  random 
phase  variations  versus  doppler.  Notice  the  3  dB  loss  now  occurs  at  about 
43  kHz  where  as  for  doppler  alone  it  occurs  at  about  60  kHz. 
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Hybrid  Array  Correlator 

A  correlator  device  which  provides  two  performance  advantage*;  over 
the  standard  digital  correlator  is  presently  beinq  developed  at  Harris  u ISP. 

The  first  performance  advantage  accrues  from  the  fact  that  the  device  operates 
on  analog  rather  than  digital  signals  so  that  the  degradation  associated  with 
one  bit  quantization  common  to  the  digital  correlator  is  avoided.  The  second 
performance  advantage  is  that  multiple  correlators,  operating  in  parallel,  are 
implemented  on  a  single  LSI  chip.  The  current  design  is  capable  of  PN  rates 
up  to  20  megachips/second,  with  growth  to  50  to  100  megachips/second  anticipated, 
and  contains  64  correlators.  A  simplified  functional  diagram  of  this  Hybrid 
Array  Correlator  (HAC)  is  shown  in  Figure  5.4-1.  As  illustrated  by  this  figure, 
the  HAC  is  a  bank  of  multiply-and-integrate  correlators  which  forms,  in  parallel, 
the  correlation  between  the  baseband  analog  signal  and  each  of  the  delayed  versions 
of  the  PN  reference.  The  delays  are  actually  achieved  by  a  64  stage  shift 
register  which  is  shifting  the  PN  reference  at  twice  the  chip  rate.  Thus, 
with  the  resulting  h  chip  offset  between  correlation  cells  the  outputs  of 
the  64  correlators  represent  the  simultaneous  search  over  32  chip  positions 
in  %  chip  steps.  The  multiplexed  output  is  therefore  the  time  sampled  pro¬ 
file,  over  32  chips  in  h  chip  steps,  of  the  correlation  between  the  input 
signal  and  the  PN  reference.  For  initial  acquisition  the  HAC  can  be  operated 
as  a  time  invariant  matched  filter  by  circulating  an  advanced  32  chip  segment 
of  the  PN  reference  through  the  shift  register.  When  the  correlation  peak 
appears  in  the  32  chip  window  the  PN  reference  timing  is  adjusted  to  place  the 
correlation  peak  in  the  center  of  the  array  and  the  shift  register  is  then 
loaded  with  the  continuous  PN  reference  sequence.  PN  code  tracking  is  accom¬ 
plished  by  simply  observing  the  location  of  the  correlation  peak  within  the 
array  and  adjusting  the  local  timing  to  keep  the  peak  centered.  Timing  need 
be  adjusted  only  often  enough  to  keep  the  peak  within  the  array,  that  is,  to 
within  ±16  chips  of  the  center. 

The  HAC  actually  contains  dual  in-phase  and  quadrature  paths  for 
both  the  input  and  PN  reference  signals,  with  each  correlator  cell  containing 
four  misers.  With 

dj(t),  dg ( t )  =  In-phase  and  quadrature  data 

I(t),  Q(t)  =  In-phase  and  quadrature  PN  sequences 

0  =  Phase  offset  between  signal  carrier  and  local 
oscillator 
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the  downconverted  I  and  Q  components  are 


Sj(t)  =  A  dj (t) I(t-T)cos  0  +  A  dg(t)Q(t-r)sin  0 
Sq( t)  =  A  dg(t)Q(t-r)cos  0  -  A  d j( t) I (t- r)si n  0 
The  four  mixers  multiply  these  signals  by  the  local  PN  sequences  to  produce 


Xj(t)  =  sr(t)  I (t) 

=  A  dj(t)I(t)I(t--r)cos  0  +  A  dg(t)I(t)Q(t-r)sin  0 
x2(t)  =  Sg(t)  Q(t) 

=  A  dg( t)Q( t)Q(t- r)cos  0  -  A  dj(t)Q(t) I (t-r)sin  0 

x3(t)  =  Sg(t)  I(t) 

=  A  dg(t)I(t)Q(t-T)cos  0  -  A  dj(t)I(t)I (t-r)sin  0 
X4(t)  =  Sj(t)  Q(t) 

=  A  djUJCKtHCt-rJcos  0  +  A  dg(t)Q(t)Q(t-T)sin  0 

These  four  signals  are  combined  and  integrated  as  appropriate  to  the  signalling 
format.  If  the  same  data  appears  on  both  the  I  and  Q  channels,  i.e. 

dj(t)  =  dq(t)  =  d(t) 
then  two  integrators  are  used  with  inputs 

Xj(t)  +  x2(t) 

and 

x4(t)  -  x3(t) 


with  the  resulting  outputs 


A  *  J  Cx1(t)  +  X2(t)]  dt 


2  A  d(t)  Rpn(t)cos  0 
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and 


=  f  [x4(t)  -  X3(t)]  dt 


2  A  d(t )  RpN(r)sin  0 


where  is  T.he  common  autocorrelation  function  of  I(t)  and  Q(t).  The 

configuration  of  each  correlation  cell  for  this  case  is  illustrated  in 
Figure  5.4-2.  If  the  data  on  the  I  and  Q  channels  are  different  then  each 
cell  can  be  configured  with  four  integrators  to  produce 


0 


~  A  dI(t)RpN(T)cos  0 
=  A  dq(t)RpN(r)cos  0 

=  A  dI(t)RpN(r)sin  0 
^  A  dQ(t)RpN(T)sin  0 


Thus,  the  HAC  approach  is  general  enough  to  permit  configuration  to  process 
BPSK,  QPSK,  OQPSK  and  MSK  with  either  one  or  two  data  streams. 

Although  the  underlying  principle  of  the  multiply-and-integrate 
correlator  has  been  known  and  applied  for  many  years,  the  design  concepts 
incorporated  in  the  HAC  permit  the  implementation  of  many  (i.e.,  64)  of  these 
correlator  cells  on  a  single  LSI  chip.  These  design  features  include  the  use 
of  on-off  switches  for  mixers,  the  use  of  shift  registers  for  obtaining  the  h 
chip  offset  between  the  cells  and  the  incorporation  of  analog  multiplexer  out¬ 
put  circuits.  The  first  two  features  are  made  possible  by  the  use  of  binary 
reference  signals.  The  use  of  binary  references  in  the  correlation  operation 
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is  appropriate  for  BPSK,  QPSK  and  OQPSK  since  the  data  and  PN  sequences  for 
these  formats  are  transmitted  as  bipolar  pulses.  For  other  spreading  forms, 
such  as  MSK,  some  correlation  loss  does  occur  by  correlating  with  a  two- level 
reference  but  generally  this  loss  is  on  the  order  of  a  few  tenths  of  a  dB. 

The  analog  output  multiplexer  feature  provides  a  convenient  means  of  reducing 
the  number  of  output  pins  to  a  practical  number.  As  presently  conceived  the 
complete  64-cell  HAC  will  be  contained  in  a  single  24-pin  IC  package. 
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SECTION  6.0 

DIRECT  SERIAL  SEARCH  VS.  MATCHED  FILTER  SEARCH 


6.0 


DIRECT  SERIAL  SEARCH  VS.  MATCHED  FILTER  SEARCH 


In  this  section  direct  serial  search  and  matched  filter  search 
techniques  are  examined  in  such  a  manner  as  to  highlight  the  procedural  and 
performance  differences  between  these  two  acquisition  processes.  DS  spreading 
is  assumed  in  order  to  minimize  implementation  complexities  which  would  tend 
to  obscure  the  key  issues. 

In  a  DS  spread  spectrum  system  the  processing  gain  is  generally 
defined  as  the  ratio  of  the  PN  chip  rate  rc  to  bit  rate  r^,  which  is  equal 
to  the  number  of  PN  chips  per  bit,  N^, 

P  =N 

PgDS  rb  Nc/b  ' 

This  processing  gain  is  usually  realized  in  the  receiver  by  a  multiply  and 
integrate  correlator  as  shown  in  Figure  6.0-1. 


Figure  6.0-1.  Multiply  and  Integrate  Correlator  Receiver 


where  the  integration  time  t-is  the  bit  length 


t  *  Tb  '  1/rb  • 

Thus  the  processing  gain  is  achieved  by  integrating  over  a  length  of  time 
which  contains  chips. 


A  similar  process  is  used  in  the  direct  serial  search  method  of  PN 
code  acquisition  illustrated  in  Figure  6.0-2. 


RECEIVED  SIGNAL 


f  <  >  * 


ACQUISITION 

DECISION 


HIT 


DISMISS 


Figure  6.0-2.  Direct  Serial  Search 


If  the  same  processing  gain  is  required  for  acquisition  as  for  data 
demodulation  then  the  product  of  the  received  signal  and  local  PN  reference  is 
integrated  for  t  seconds  (Nc^  chips)  and  the  result  used  for  making  an  acqui¬ 
sition  decision,  usually  by  a  threshold  device.  If  a  "hit"  occurs,  a  validation 
procedure  is  initiated;  if  a  "no  hit"  decision  occurs  the  PN  reference  is  slewed 
^-chip  and  the  process  repeated.  The  basic  limitation  to  search  speed  is  that 
the  correlator  integrates  over  a  completely  new  set  of  Nc^b  chips  of  the  received 
signal  for  each  test.  Thus  the  search  rate  is 

Search  Rate  for  r 

Direct  =  =  r.,  %-chip  positions  per  second. 

Serial  Search  Nc/b  b 

The  search  rate  can  be  significantly  increased  by  replacing  the 
multiply  and  integrate  correlator  by  an  MF  device.  Acquisition  of  DS  by  an  MF 
is  best  illustrated  by  considering  a  digital  correlator  of  Figure  6.0-3.  The 
received  signal  is  sampled  and  digitized  (to  one  bit)  at  twice  the  chip  rate 
and  serially  loaded  into  a  shift  register.  An  advance  segment  of  the  PN  refer¬ 
ence  is  loaded  into  a  holding  register.  The  contents  of  the  two  registers  are 
correlated  by  comparing  the  contents  of  the  corresponding  stages,  generating  a 
"+1"  if  the  contents  match  and  a  "-1"  if  they  do  not,  and  summing  the  results. 

The  correlation  value  is  tested  at  the  clock  rate  (which  is  twice  the  PN  chip 
rate)  so  that  the  search  rate  for  this  method  is 

Search  Rate  for  _  2  ij-chip  positions  per  second. 


CORRELATION 


Figure  6.0-3.  MF  For  DS  Acquisition 


This  indicates  a  reduction  in  search  time  over  the  direct  serial  search  method 
by 


Matched  Filter  Search  Rate 


Direct  Serial  Search  Rate 


=  2  N 


c/b 


Therefore  if  a  processing  gain  of  30  dB  is  required  during  acquisition  then 
the  MF  reduces  the  time  to  search  over  the  total  uncertainty  region  by  a 
factor  of  2000.  The  improvement  in  search  time  is  a  consequence  of  the  fact 
that  only  a  new  ^-chip  of  received  signal  is  required  for  each  correlation  test 
since  prior  samples  are  already  stored  in  the  shift  register.  (As  a  practical 
note,  the  processing  gain  of  an  MF  is  proportional  to  its  length  and  at  present 
no  digital  correlator  is  available  which  will  provide  30  dB  of  processing  gain 
in  a  single  package.) 

A  basic  different  in  receiver  operations  for  the  multiply  and  inte¬ 
grate  and  MF  methods  is  that  in  the  multiply  and  integrate  receiver  the  local 
PN  generator  is  running  continuously,  whereas  in  the  MF  receiver  the  PN  segment 
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in  the  holding  register  is  fixed.  The  multiply  and  integrate  acquisition 
process  may  be  visualized  as  having  both  the  received  and  local  PN  sequences 
moving  together  at  the  chip  rate,  with  ’2-chip  shifts  in  the  relative  phase 
introduced  every  t  seconds.  These  ^-chip  shifts  are  continued  until  alignment 
of  the  two  sequences  occurs.  In  the  MF  receiver  the  local  PN  reference  is 
stationary  and  the  received  PN  sequence  slides  past  until  alignment  occurs, 
at  which  time  the  local  PN  generator  is  enabled.  These  processes  are  illus¬ 
trated  in  Figure  6.0-4. 
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Figure  6.0-4(a).  Direct  Serial  Search 
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Figure  6.0-4(b).  Matched  Filter  Search 
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7.0  COMPARISON  OF  OS,  FH  AND  FH/DS 

In  Section  6.0  acquisition  performance  was  examined  as  a  function 
of  the  structure  of  the  receiver's  acquisition  system,  i.e.,  direct  search 
versus  MF.  In  this  section  acquisition  performance  is  examined  as  a  function 
of  the  type  of  spreading  function  employed. 

7. 1  Direct  Sequence 

The  results  of  Section  6.0  can  be  used  directly  in  characterizing 
DS  acquisition  performance  relative  to  the  other  forms.  For  the  direct  search 
method  the  search  rate  is  the  data  rate. 


Search  Rate  for  ,  .  . 

Direct  Search  =  rb’  "2-chlp  P°s1t1on  Per  seco"d 
whereas  for  the  MF  method  the  search  rate  is  twice  the  PN  chip  rate. 


Search  Rate 
for  MF 


=  2  rc,  *s-chip  position  per  second. 


In  both  methods  the  search  is  normally  conducted  in  ?s-chip  increments  so  that 
when  synchronization  is  achieved  timing  uncertainty  between  the  transmitter 
and  receiver  is  resolved  to  within  Jj-chip. 

7.2  Frequency  Hop 

In  frequency  hop  (FH)  systems  processing  gain  is  a  function  of  the 
number  of  hop  frequencies  rather  than  the  number  of  PN  chips  per  bit.  With 


=  Bandwidth  of  the  FH  spread  spectrum  signal 
B<-  =  Bandwidth  of  each  frequency  slot 
Np  =  Number  of  hop  frequencies 
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then  the  processing  gain  is 


^SS 

pgFH  ”b^~  =  nf 


The  basic  FH  transmitter  is  as  shown  in  Figure  7.2-1. 


Figure  7.2-1.  Basic  FH  Transmitter 


The  PN  generator  is  used  to  generate  the  binary  address  of  the  hop  frequency 
by  providing  one  new  chip  to  the  shift  register  for  each  hop.  If  there  are 
NH/b  fre(Iuency  hops  Per  bit>  then  the  rate  the  generator,  >"cpH.  is 
related  to  the  bit  rate  by 

rcFH  =  NH/b  rb 

The  corresponding  rate  of  the  PN  generator  in  a  DS  system  is 
rcDS  =  Nc/b  rb 
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Since  there  are  typically  many  orders  of  magnitude  more  chips  per  bit  in  a  DS 
system  than  there  are  frequency  hops  per  bit  in  a  FH  system,  the  ratio  of  the 
two  PN  rates, 

rcDS  =  Nc/b 
rcFH  NH/b 

is  typically  many  orders  of  magnitude.  Thus  the  width  of  a  FH  chip,  TcpH,  is 
many  orders  of  magnitude  greater  than  the  width  of  a  DS  chip,  T 

7.2.1  Direct  Search  of  FH 

Direct  search  for  FH  is  conducted  in  a  manner  analogous  to  DS  direct 
search,  as  illustrated  in  Figure  7. 2. 1-1. 


Figure  7.2. 1-1.  Direct  Search  for  FH  Acquisition 

By  searching  at  a  rate  equal  to  the  bit  rate,  r^,  and  slewing  the  PN  generator 
in  Vchip  steps  this  direct  search  progresses  through  the  timing  uncertainty 
region  at  the  same  basic  rate  and  in  the  same  manner  as  the  DS  direct  search. 
But  with  the  greater  width  of  the  FH  chip  a  larger  portion  of  the  uncertainty 

region  is  covered  in  each  step,  by  the  factor 
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1  T  M 

2  cFH  _  c/b 

IT  N 

2  cDS  NH/b 


Therefore  the  time  required  to  search  through  a  given  uncertainty  range  is 

reduced  by  the  factor  ^  course  we  canlt  e*Pect  something  for 

nothing  —  the  price  is  that  when  timing  is  resolved  to  within  ^s-chip  in 

each  of  the  two  systems  the  resolution  is  coarser  in  the  FH  system  by  the 

same  N  factor.  This  would  seem  to  indicate  that  FH  is  an  attractive 

c/b  H/b 

acquisition  aid,  rapidly  reducing  a  large  initial  timing  uncertainty  to  a 
range  which  permits  rapid  acquisition  of  DS  by  a  MF.  However,  it  must  be  noted 
that  when  the  frequency  hop  rate  is  one-half  the  bit  rate,  both  FH  acquisition 
and  MF  acquisition  of  DS  show  the  same  acquisition  speed  advantage  over  direct 
search  of  DS,  and  the  MF  acquisition  of  DS  results  in  finer  timing  resolution. 
For  lower  FH  rates  FH  does  show  a  speed  advantage  over  the  MF  acquisition  of 
DS  but  an  even  greater  reduction  in  resolution. 


7.2.2  MF  for  FH 

The  obvious  next  step  is  to  reduce  the  time  required  for  acquisition 
of  FH  by  utilizing  a  FH  MF  such  as  in  Figure  7. 2. 2-1. 


Figure  7. 2. 2-1.  MF  for  FH  Acquisition 


7-4 


|  A'V’W*** 


I 

f 


In  this  receiver  the  PN  generator  is  held  at  an  advanced  state  so  that  the 
Discrete  Frequency  Synthesizer  produces  a  constant  frequency.  The  input  to 
the  N  parallel  channels  is  therefore  a  FH  pattern  which  is  simply  offset  in 
frequency  from  the  received  FH  pattern.  The  output  of  a  particular  BPF  is 
large  when  the  frequency  of  the  input  pattern  matches  the  frequency  of  the 
filter.  The  filter  frequencies  and  the  corresponding  delays  are  selected  to 
match  the  FH  pattern  which  will  be  generated  by  the  PN  sequence  prior  to  and 
including  the  PN  segment  held  in  the  PN  generator.  A  typical  advanced  FH 
pattern  is  shown  in  Figure  7. 2. 2-2. 
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Figure  7. 2. 2-2.  Advanced  FH  Pattern 
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With  this  arrangement  the  time-hopped  pulses  at  the  envelope  detector  outputs 
will  be  in  alignment  at  the  input  to  the  summing  network,  providing  a  S/N  ad¬ 
vantage  over  and  above  the  FH  processing  gain.  Thus  if  the  S/N  at  the  acqui¬ 
sition  detector  in  this  system  is  to  be  maintained  at  the  same  level  as  in  the 
conventional  FH  direct  search  system,  the  S/N  ratio  of  each  channel  can  be 
reduced.  This  permits  the  BPF  bandwidths  to  be  increased  and  the  sampling  rate 
to  be  increased.  Ideally  (ignoring  non-coherent  comining  loss)  these  parameters 
are  linearly  related,  that  is,  the  S/N  ratio  into  the  acquisition  detector  is 
unchanged  if  the  number  of  channels  is  increased  from  1  to  N,  the  bandwidths 
increased  by  N  and  the  sampling  rate  increased  by  N.  In  general  the  number  of 
parallel  channels,  N,  will  have  to  be  somewhat  larger  than  the  number  of  samples 


I 
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per  frequency  hop  interval,  N<~,  to  overcome  the  non-coherent  combining  loss. 

With  samples  per  FH  interval,  the  undelayed  samples  in  each  channel  will 
be  as  shown  in  Figure  7. 2. 2-3  and  the  delayed  samples  as  in  Figure  7. 2. 2-4. 

The  acquisition  detector  will  be  triggered  at  the  time  of  first  coincidence  of 
the  signal  which  occurs  no  later  than  T^/N^  seconds  after  the  beginning  of  the 
Nth  interval.  Thus  acquisition  with  this  MF  approach  provides  timing  resolution 


Resolution  of  MF 
Acquisition  of  FH 


which  for  large  Ns  would  indicate  a  considerable  improvement  over  the  TcFH/2 
resolution  of  direct  search  acquisition  of  FH.  In  fact,  by  noting  that 


cFH  = 
N- 


we  see  that  MF  acquisition  of  FH  can  achieve  the  TcDS/2  resolution  of  DS  by 
sampling  at  a  rate  of 

Ns  =  j  N£/fr 

NH/b 

samples  per  frequency  hop.  This  is  a  surprising  result  since  the  MF  approach 
was  anticipated  to  produce  an  improvement  in  acquisition  time,  not  timing  re¬ 
solution. 

The  question  now  is:  Does  the  MF  approach  to  FH  acquisition  provide 
any  improvement  in  acquisition  time?  To  answer  this  question  recall  that  it 
was  noted  earlier  that  direct  search  of  FH  results  in  the  same  acquisition 
time  performance  as  MF  acquisition  of  DS  when  the  frequency  hopping  rate  is 
one-half  the  bit  rate,  and  FH  gains  some  advantage  when  the  hop  rate  is 
reduced.  Thus  the  question  can  be  answered  by  comparing  the  search  rates  of 
MF  acquisition  of  FH  and  MF  acquisition  of  DS: 


Search  Rate  for  MF 
Acquisition  of  FH 


1  FH-chip  of  time  searched 
per  FH  time  interval. 

7-6 


t -  ,  FH-chips  of  time  searched  per  second 

cFH 


^c/b^H/b  ,  DS-chips  of  time  searched  per 
TcFH  second 


^c/b^H/b  ,  DS-chips  of  time  searched 


WVb^cDS  Per  second 


_1 _  ,  DS-chips  of  time  searched  per  second 


r  ,  DS-chips  of  time  searched  per  second 


Search  Rate  for  MF 
Acquisition  of  DS 


2  r  n<.  ,  \  DS-chips  of  time  searched  per 
cu*  second 


*  rcDS  ’  D^-chips  of  time  searched  per  second 

Thus  for  any  initial  timing  uncertainty  the  two  MF  approaches  would  achieve  acqui¬ 
sition  in  the  same  time.  In  addition,  if  the  sampling  rate  is  high  enough  in  the 
FH  MF  system  the  timing  resolutions  at  acquisition  are  the  same  for  the  two  sys¬ 
tems.  The  relationship  between  the  four  acquisition  methods  can  be  summarized  by: 


REDUCED  ACQ. 
TIME,  EQUAL 
RESOLUTION 


ME  SEARCH 
OF  DS 


REDUCED  ACQ.  TIME, 
REDUCED  RESOLUTION 


EQUAL  ACQ.  TIME, 

IMPROVED 

RESOLUTION 


MF  SEARCH 
OF  FH 


7.3 


In  FH/DS  systems  the  total  processing  gain  (in  dB)  is  approximately 
equal  to  the  sum  of  the  processing  gains  due  to  the  FH  spreading  and  the  DS 
spreading. 


pgFH/DS^dB^  "  pgFH^dB^  +  pgDS^dB^  ’ 

In  terms  of  the  nomenclature  defined  on  the  basic  FH/DS  transmitter  block 
diagram  of  Figure  7.3-1,  the  processing  gain  is 

pgFH/DS  =  10  log  NF  +  10  log  Nc/b  •  (7.3-1) 

This  processing  gain  is  essentially  fixed  by  the  bandwidth  of  the  spread 
spectrum  signal  and  may  be  divided  between  the  two  spreading  functions  as 
desired.  The  frequency  hopping  rate  (number  of  frequency  hops  per  bit,  N^b) 
does  not  enter  into  the  processing  gain  per  se  but  rather  is  determined  by 
considering  the  frequency-following  capability  of  the  jamming  threat.  Typi¬ 
cally  the  number  of  DS  chips  per  bit,  Nc^,  is  two  or  more  orders  of  magnitude 
greater  than  the  number  of  frequency  hops  per  bit,  N^b,  so  that  the  DS  PN  rate 
is  correspondingly  greater  than  the  FH  PN  rate. 

7.3.1  Direct  Search 

The  receiver  structure  for  direct  search  acquisition  of  FH/DS  is 
shown  in  Figure  7. 3. 1-1.  It  is  assumed,  as  before,  that  the  multiply  and 
integrate  correlator  for  the  DS  despreading  integrates  over  a  bit  interval  so 
that  the  search  rate  is  equal  to  the  bit  rate.  The  PN  generator  is  stepped 
through  the  timing  uncertainty  region  in  increments  equal  to  h  of  the  smaller 
chip  width,  i.e.,  either  4  DS  chip  or  \  FH  chip  depending  on  which  is  smaller. 
As  noted  earlier  the  DS  PN  rate  is  usually  much  greater  than  the  FH  PN  rate 
so  that  the  DS  chip  width  is  usually  smaller  and  the  timing  uncertainty 
region  is  searched  in  increments  of  %  a  DS  chip.  Consequently,  once  acquisi¬ 
tion  is  achieved  the  timing  uncertainty  is  resolved  to  within  H  a  DS  chip. 

Note  that  for  systems  with  the  same  processing  gain  pure  FH  typically  has  the 
lowest  chip  rate,  pure  DS  has  the  highest  and  FH/DS  is  in-between.  Therefore, 
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FH  has  the  largest  search  increment,  fastest  acquisition  and  coarsest  timing 
resolution;  DS  has  the  smallest  search  increment,  slowest  acquisition  and 
finest  timing  resolution;  and  FH/DS  is  in-between  in  all  three  respects. 

Thus  FH/DS  provides  a  means  of  compromising  between  acquisition  time  and 
timing  resolution. 

More  quantitatively  let  AT  be  the  timing  uncertainty,  assume  the 
search  is  conducted  in  Vchips  increments  at  a  rate  equal  to  the  data  rate, 

Rb,  and  let  Rc  denote  the  chip  rate  appropriate  to  the  type  of  spreading 
being  considered.  (There  is  no  ambiguity  regarding  Rc  in  pure  DS  or  pure 
FH  systems,  in  FH/DS  systems  the  appropriate  Rc  is  the  higher  of  the  DS  and 
FH  chip  rates.)  Then  the  total  number  of  Jj-chip  increments  in  the  uncertainty 
region  is 

#*s-chip 

Increments  of  =  2  AT  R 
Uncertainty  c 

and  the  acquisition  time  (or,  more  precisely,  the  time  for  one  complete  search) 
is 

r\  .*r 

Acquisition  Time  =  6  R 

Kb  c 

The  timing  uncertainty  remaining  at  acquisition  is  %-chip. 

Timing  Uncertainty  _  v  T  v  _1_ 
at  Acquisition  2  c  ^  Rc  ’ 

Typically  R  increases  in  going  from  FH  to  FH/DS  to  DS.  In  DS  systems  R  is 
c  c 

set  by  the  required  processing  gain;  in  FH  systems  by  the  frequency  hopping 
rate  required  to  avoid  a  frequency-following  jammer  (with  processing  gain 
coming  from  the  number  of  hop  frequencies).  FH/DS  systems  provide  an  addi¬ 
tional  degree  of  freedom  in  that  the  three  parameters,  DS  chip  rate,  FH  rate 
and  number  of  hop  frequencies,  can  be  selected  as  a  set  to  simultaneously  pro¬ 
vide  the  required  AJ  performance  and  reduced  acquisition  time.  The  performance 
comparisons  for  direct  search  can  be  summarized  by: 


Spreading 

Search 

Increment 

Size 

Acquisition 

Time 

Timing 

Resolution 

FH 

1 

1 

1 

FH/DS 

Decreasing 

Increasing 

Increasing 

DS 

L 

i 

1 

7.3.2  Matched  Filter 

In  the  MF  method  of  acquiring  FH/DS,  as  in  the  MF  method  for  all 
types  of  spreading  functions,  the  PN  generator  is  fixed  at  some  advanced  state 
and  the  receiver  simply  awaits  coincidence  between  the  received  and  local 
spreading  functions.  The  receiver  structure,  as  shown  in  Figure  7. 3. 2-1,  is 
the  same  as  for  direct  search  of  FH/DS  except  that  the  integrate  and  multiple 
correlator  for  the  DS  despreading  is  replaced  by  a  MF. 

Comparing  MF  acquisition  of  FH/DS  with  MF  acquisition  of  DS,  the 
FH/DS  system  has: 

•  Lower  DS  Chip  Rate 

•  Larger  Search  Increments 

(i.e.  Larger  Jj-Chip  Width) 

•  Lower  Search  Rate 

(i.e.  Twice  the  DS  Chip  Rate) 

The  larger  search  increments  dictate  coarser  timing  resolution  at  acquisition. 
Since  the  increase  in  size  of  the  search  increment  is  exactly  offset  by  the 
decrease  in  search  rate  the  total  search  time,  or  acquisition  time,  remains 
unchanged  from  the  DS  case.  Comparison  of  acquisition  performances  for  the 
six  possible  combinations  of  spreading  function  and  receiver  techniques  is 
summarized  in  Figure  7. 3. 2-2. 

Intuitively  it  would  seem  that  a  parallel  bank  of  channels  at 
different  frequencies  (analogous  to  the  FH  system  of  Figure  7. 2. 2-1,  not  parallel 
search  operations)  would  improve  the  timing  resolution  in  the  FH/DS  system  as 
it  did  in  the  FH  system.  However,  the  resolution  in  the  FH/DS  is  limited  by 
the  fact  that  the  received  signal  is  sampled  and  correlated  with  the  reference 
DS  signal  at  twice  the  DS  chip  rate.  Therefore  timing  resolution  can  be  no 
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Summary  of  Acquisition  Performances 


better  than  h  a  DS  chip.  Likewise  multiple  channels  can  not  improve  the 
acquisition  time  since  the  search  rate  is  set  by  the  rate  of  the  incoming 
signal.  Multiple  channels  would  provide  a  higher  S/J  at  the  input  to  the 
acquisition  detector  which  would  improve  and/or  Pp^,  parameters  which 
have  not  been  considered  in  the  performance  comparisons  presented  in  this 
section. 


SECTION  8.0 

BASIC  ACQUISITION  STRATEGIES 


8.0  BASIC  ACQUISITION  STRATEGIES 

In  this  part  we  turn  our  attention  to  the  basic  strategies  for 
transferring  timing  information  from  transmitter  to  receiver.  The  context 
is  military  spread  spectrum  applications  which  imposes  the  additional  require¬ 
ment  that  protection  against  intelligent  jamming  not  be  unduly  compromised 
during  the  acquisition  process.  The  intent  is  to  characterize  the  basic 
acquisition  strategies  available,  rather  than  specific  acquisition  techniques 
or  implementations,  in  such  a  manner  as  to  permit  order-of -magnitude  compari¬ 
sons  in  potential  acquisition  time  and  to  assess  to  some  degree  the  potential 
vulnerability  to  intelligent  jamming. 

The  "acquisition  time"  expressions  developed  are  actually  expressions 
for  "time  to  complete  one  search  of  the  timing  uncertainty  region"  or  "time  to 
transmit  timing  information,"  etc.  rather  than  more  precise  characterizations 
of  acquisition  time  such  as  "mean  time  to  acquire."  Effects  of  false  alarms 
and  missed  detections  have  been  neglected  so  that  the  key  parameters  which 
dictate  performance  potential  can  readily  be  identified.  Likewise,  timing 
uncertainty  is  expressed  in  terms  of  the  number  of  PN  chips,  L,  of  uncertainty 
which  exists  between  the  transmitter  and  receiver  in  order  to  minimize  the  num¬ 
ber  of  parameters.  L  may  be  the  total  length  of  the  PN  code  or  some  portion  of 
it  depending  on  the  extent  of  timing  uncertainty.  Timing  uncertainty  in  seconds, 
AT,  is  related  to  L  by  the  PN  chip  rate,  Rc, 

L  =  AT  Rc  (8.0-1) 

For  a  given  value  of  AT  the  corresponding  value  of  L  can  vary  greatly  depending 
on  the  type  of  spreading  employed;  for  example  for  the  same  AT,  L  would  typically 
be  several  orders  of  magnitude  greater  for  a  direct  sequence  system  than  for  a 
frequency  hop  system  due  to  the  typically  large  difference  in  PN  rates.  These 
are  systems  considerations  which  are  implicitly  acounted  for  in  the  acquisition 
models  presented,  but  which  are  not  explicitly  shown  in  order  that  the  key  fea¬ 
tures  of  the  basic  strategies,  not  fundamental  system  concepts,  be  emphasized. 

Basically  with  a  PN  code  of  length  L  chips  timing  information  can  be 
transferred  from  transmitter  to  receiver  by  either  transmitting  the  PN  code- 
driven  spreading  function  and  searching  for  the  correct  code  phase  in  the 
receiver  or  by  transmitting  the  code  phase  as  an  N-bit  binary  number  using  a 
fixed  spreading  function,  where 
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The  first  method  has  the  advantage  of  full  AJ  protection  during  acquisition 
but  results  in  excessively  long  acquisition  times  for  typically  long  code 
lengths.  The  second  method  results  in  very  rapid  acquisition  but  is  quite 
vulnerable  to  smart  jammers  since  the  spreading  function  remains  fixed  while 
transmitting  the  N  bits  of  timing  information.  A  third  approach  is  to  use 
a  combination  of  these  two,  that  is  partially  resolving  the  timing  uncertainty 
by  transmitting  and  searching  over  a  short  PN  code-driven  spreading  function 
and  then  transmitting  the  remaining  timing  information  as  data  on  the  short 
code-spread  signal. 

8. 1  Search  Methods 

There  are  two  basic  search  methods  of  acquisition,  direct  search 
(also  commonly  referred  to  as  conventional  serial  search)  and  matched  filter 
(MF)  search.  In  direct  search  the  receiver  correlates  the  incoming  signal 
with  the  local  reference  over  a  t  second  time  interval  and  makes  a  decision 
as  to  whether  or  not  the  two  are  synchronized.  If  not  then  the  timing  of  the 
local  reference  is  incremented  one-half  chip  and  tested  again.  This  continues 
until  synchronization  is  achieved.  If  there  are  S  such  search  devices  each 
independently  searching  1/S  of  the  uncertainty  region,  then  the  time  required 
to  search  the  entire  2L  ^-chip  element  uncertainty  region  is 


The  search  rate  for  this  technique  is 

DS:  R  =  1/t,  *5-chip  positions/second  so 
s  that  Tg,.  may  be  expressed  as 

TDS  =  lrs  (8.1-1) 

In  the  MF  search  method  the  reference  section  of  an  MF  is  loaded 
with  an  advanced  segment  of  the  incoming  spreading  waveform,  checks  the  corre¬ 
lation  between  the  incoming  signal  and  the  frozen  reference  at  twice  the  chip 
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rate  and  signals  synchronization  when  the  correlation  exceeds  a  threshold 
value.  The  search  rate  in  this  case  is  thus 


MF:  =  2  Rc,  ^-chip  position/second 


With  S  MF's  each  MF  is  assigned  (2L/s)  %-chip  positions.  The  time  required  to 
send  over  the  complete  uncertainty  region  is  the  time  required  to  initially  load 
the  devices  and  obtain  the  first  correlation  test  values,  t  ,  plus  the  time  for 
each  MF  to  search  over  the  remaining  (2L/S)  positions  at  a  rate  of  R$,  i.e.. 


TMFS 


(T. 1-2) 


The  initial  response  time,  tq,  is  a  function  of  the  length  of  the  MF, 


^MF  £  ^F  length  in  chips 


and  the  i 


nput  chip  rate  R  , 
c 


T 


o 


=  2Aif 


(8.1-3) 


Note  that  the  two  expressions  for  search  time  are  similar  even  though 
the  actual  search  times  may  be  quite  different.  The  search  rate  for  the  MF  is 
twice  the  chip  rate;  for  direct  search  it  is  usually  on  the  order  of  the  bit 
rate  which  is  typically  orders  of  magnitude  lower  than  the  chip  rate.  There¬ 
fore  the  MF  usually  provides  orders  of  magnitude  faster  acquisition  than  does 
direct  search.  Nonetheless  a  common  expression  can  be  obtained  for  T^  and  T^p<- 
by  manipulating  Tnc  as  follows: 


But  t  is  the  time  required  for  each  correlation  test,  including  the  first  one 
so  that  for  direct  search 


and 


To  T 


T 


0 


which  is  the  same  as  the  MF  expression.  Therefore,  a  common  expression  for 
the  search  methods  is 


(8.1-4) 


8.2  The  Data  Transmission  Method 

An  alternative  to  the  search  method  of  acquisition  is  the  trans¬ 
mission  of  timing  information  as  N  bits  of  data,  where 

L  <  2N 

using  a  fixed  spreading  function.  The  spreading  function  may  be  a  fixed  segment 
of  a  PN  sequence,  a  fixed  chirp  waveform,  etc.  Transmitting  the  N  bits  of  timing 
information  at  a  bit  rate  of  bits  per  second  results  in  an  acquisition  time 
(i.e.,  time  to  transmit  the  N  bits)  of 

Tqt  =  ■—  ,  seconds.  (8.2-1) 

The  ratio  of  acquisition  times  for  the  data  transmission  method  and  the  search 
method  is 


%  SN 

rd  K 
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which  for  large  N  indicates  many  orders  of  magnitude  faster  acquisition  by 
the  data  transmission  method,  Thisimproved  performance  is  of  course  due  to 
the  fact  that  there  is  no  uncertainty  in  the  spreading  function,  which  also 
makes  the  receiver  quite  vulnerable  to  smart  jamming. 

8.3  The  Combined  Search  and  Data  Transmission  Method 

Basically,  the  move  from  a  search  method  of  acquisition  time  to  the 
data  transmission  has  been  one  from  great  uncertainty  in  the  spreading  function, 
with  correspondingly  large  acquisition  time  but  considerable  immunity  to  smart 
jamming,  to  no  uncertainty  in  the  spreading  function,  with  correspondingly  small 
acquisition  time  but  considerable  vulnerability  to  small  jamming.  A  combination 
of  these  two  methods  provides  a  reasonable  compromise  between  these  two  extremes, 
and,  under  certain  conditions,  results  in  an  even  somewhat  faster  acquisition 
than  even  the  data  transmission  method.  Many  so  called  rapid  acquisition  tech¬ 
niques  employ  some  variation  of  this  idea. 

In  this  combined  method  acquisition  of  a  long  PN  code  with 


LLC  2 


chips — bits  of  timing  information 


is  achieved  by:  (1)  transmitting  a  short  PN  code  -  driven  spreading  function 
with 


LSC  ~  2  chips— bits  of  timing  information 


(2)  acquiring  this  short  code  by  a  search  method  and  (3)  transmitting  the  remain¬ 
ing  (nlc~NS(0  bits  of  timing  information  as  data  on  the  short  code  spread  signal. 
The  time  to  acquire  using  this  method  therefore  is  the  time  to  search  over  Lsc 
chips  plus  the  time  to  transmit  bits  of  data.  The  time  to  conduct  one 

search  over  the  short  code  is 


T$  =  l  .  j)  „  T( 


and  the  time  to  transmit  ( )  bits  of  data  is 


so  that  the  total  time  to  acquire  the  long  PN  code  using  this  combined  method 
is 


(8.3-1) 


The  performance  of  this  method  relative  to  the  other  two  is  best  demonstrated 
by  an  example. 
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Example 

Consider  the  acquisition  of  a  long  code-driven  DS  system  with 
N  =  50 

L  =  1.13  x  1025  chips 
Rc  =  10  megachips/second 
S  =  10 


Direct  Search 

Assume  that  the  nominal  data  rate  for  the  system  is  5  Kbps  and  that 
acquisition  must  be  accomplished  with  the  same  processing  gain.  Then  the  search 
rate  is 


3 

R$  =  5  x  10  ,  *s-chip  positions/second 
and  the  time  for  one  search  over  the  entire  long  code  is 
Tj  =  4.5  x  10*°  seconds. 
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Obviously,  the  direct  search  method  is  unsuitable  for  initial  acquisition  of 
the  long  code. 

MF  Search 

With  an  MF,  the  search  rate  is  twice  the  chip  rate  so  that 
Rg  =  20  x  10^,  J^-chip  positions/second 


and 


T<.  =  11.3  x  10^  seconds. 

The  MF  thus  provides  a  considerable  improvement  (approximately  4000  times)  in 
acquisition  time  but  is  still  unsuitable  for  this  initial  acquisition  situation. 

Data  Transmission 

With  the  50  bits  of  timing  information  transmitted  at  a  data  rate 
of  5  Kbps,  the  data  transmission  method  requires  only 

Tdt  =  10  milliseconds 

This  is  a  tremendous  improvement  in  acquisition  time  but  with  considerable  vul¬ 
nerability  to  smart  jamming  since  the  spreading  function  is  fixed. 


Combined  Method 

The  time  to  acquire  using  the  combined  method  is  a  function  of  the 
short  code  length. 


T 


C 


1 

20xl06 


(50-Nsc) 

5xl03 


(8.4-1) 


which  is  plotted  in  Figure  8.4-1.  It  is  noted  that  a  minimum  occurs  and  the 
value  of  7.4  milliseconds  (corresponding  to  a  short  code  length  of  approximately 
28,000  chips)  actually  represents  an  improvement  over  the  data  transmission 
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method.  This  raises  some  interesting  questions  regarding  optimization  of 
the  combined  method  and  its  general  performance  relative  to  that  of  the  data 
transmission  method. 


8.5  The  Optimum  Combined  Method 

The  optimum  value  of  the  short  code  length  in  the  combined  method 
can  be  determined  by  expressing  as 

N$C  =  log2  LSC 


in  the  expression  for  Tc>  considering 
ferentiating  Tc  with  respect  to  l_sc. 


Lsc  as  a  continuous  variable  and  dif- 
The  resulting  optimum  value  is 


L  -  5 

LSC  OPT  ( 21n2 )  RQ 


(8.5-1) 


and 


S  Rs 

NSC  OPT  log2  (21n2)  Rd 


(It  is  interesting  to  note  that  the  optimum  short  code  length  is  not 
a  function  of  the  long  code  length.)  The  corresponding  optimum  acquisition 
time  using  the  combined  method  is 


*  Ro  (v  *l) 

+  \c  -  1o92  jVsinfy  ftD  j  | 

For  the  example  considered  earlier 

For  the  example  considered  earlier 

For  the  example  considered  earlier 
L$c  opt  =  28,854  chips 

Tc  Qpj  =7.4  milliseconds 

which  match  the  values  indicated  in  Figure  8.4-1. 
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(8.5-2) 
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I 

|  Figure  8.4-1.  Performance  of  the  Combined  Method 
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The  next  question  of  interest  is:  Under  what  conditions  will  the 
combined  method  be  faster  than  the  data  transmission  method?  Since  (1)  T^. 
approaches  the  TDT  value  as  LSC~~0,  (2)  approaches  the  large  T<.  value 
when  LS£  —  L^,  and  (3)  vs.  has  a  unique  minimum,  then  can  be  less 
than  TDT  only  if  the  slope  of  T^.  vs.  is  negative  at  the  origin,  i.e.. 


dT 


dL 


SC 


LSC  0 


This  leads  to  the  condition 


SRs  >  (2  In  2)  Rq  =  1.39  RQ. 


(8.5-3) 


Note  that  R$  is  the  search  rate  of  each  of  the  S  independent  search  devices  so 
that  SR$  is  the  total  search  rate.  The  requirement  therefore  is  that  the  total 
search  rate  for  the  short  code  search  be  somewhat  larger  than  the  data  rate  for 
the  data  transmission  method.  R$  is  usually  several  orders  of  magnitude  greater 
than  TD  for  an  MF  search  implementation  so  that  the  combined  method  would  usually 
be  faster  than  the  data  transmission  method  if  an  MF  is  employed.  And,  of 
course,  the  main  advantage  of  the  combined  method  is  the  increased  protection 
against  smart  jamming. 

8.6  Conclusions 

As  is  widely  recognized,  an  MF  can  provide  a  significant  reduction 
in  spread  spectrum  signal  acquisition  time  over  direct  search  techniques  and 
for  reacquisition  or  acquisition  with  small  timing  uncertainty  may  result  in 
acceptable  acquisition  performance.  However,  for  acquisition  with  large  and 
possibly  even  total  timing  uncertainty  even  the  MF  performance  may  be  accep¬ 
table  by  many  orders  of  magnitude.  In  such  cases  the  only  alternative  is  to 
reduce  the  randomness  in  the  spreading  function  and  tradeoff  smart  jammer 
protection  for  decreased  acquisition  time.  This  basic  conflict  arises  from 
the  fact  that  any  randomness  placed  in  the  spreading  function  to  confuse  a 
smart  jammer  also  confuses  an  unsynchronized  receiver.  Completely  removing 
this  randomness,  as  in  the  data  transmission  method,  results  in  very  rapid 


acquisition  but  also  makes  the  system  very  vulnerable  to  smart  jamming. 
The  combined  method  presents  a  means  for  compromising  between  the  two. 
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CODE  VALIDITY  INTERVAL  TECHNIQUES 

The  acquisition  problem  can  obviously  be  circumvented  by  using 
fixed  spreading  waveforms;  for  example,  by  using  a  different  but  fixed 
segment  of  a  PN  sequence  to  represent  each  symbol  in  the  data  symbol  alphabet. 
The  receiver  is  then  simply  a  bank  of  time  invariant  matched  filters,  one 
matched  to  each  symbol  sequence,  and  relative  timing  between  transmitter  and 
receiver  becomes  completely  irrelevant  insofar  as  despreading  is  concerned. 

But  obviously  the  use  of  fixed  spreading  waveforms  makes  the  system  quite 
vulnerable  to  jamming  by  intelligent  jammers.  This  vulnerability  can  be 
reduced  in  systems  with  limited  timing  uncertainty  by  periodically  changing 
the  spreading  functions.  The  trick  is  to  change  the  set  of  spreading  functions 
at  a  rate  which  is  fast  enough  to  confuse  the  jammer  but  slow  enough  so  that 
both  the  transmitter  and  receiver  are  using  the  same  set  of  spreading  functions 
"most  of  the  time."  Thus  the  trade-off  is  between  speed  of  the  smart  jammer 
and- timing  accuracy  in  the  communications  system.  The  time  interval  over 
which  each  set  of  spreading  waveforms  is  valid  is  referred  to  as  the  Code 
Validity  Interval  (CVI). 

9.1  Hard  Edge  CVI 

In  conventional  (so-called  hard  edge)  CVI  systems  the  waveforms 
are  fixed  throughout  each  CVI  and  changed  abruptly  at  the  end  of  each  CVI. 

Thus  the  MFs  in  the  receiver  are  exactly  matched  to  the  transmitter’s  set  of 
waveforms  during  that  part  of  the  CVI  in  which  both  the  transmitter  and  receiver 
are  using  the  same  set,  and  not  matched  during  that  part  in  which  they  are  using 
adjacent  sets  due  to  timing  errors.  The  problem  occurs  at  the  beginning  and 
end  of  each  CVI  so  that  data  loss  can  be  avoided  by  providing  a  guard  band 
around  the  transition  point  during  which  no  transmissions  occur.  However, 
inclusion  of  a  guard  band  reduces  the  available  transmission  time  so  that  the 
transmitted  symbol  rate  must  increase  (if  the  data  rate  is  to  remain  constant) 
and  the  processing  gain  decreases  correspondingly.  For  example,  to  protect 
against  data  loss  in  a  system  with  maximum  timing  offset  between  transmitter 
and  receiver  of  ±30%  of  the  CVI  would  require  a  total  guard  band  of  60%  of 
the  CVI.  This  would  require  a  2. 5  times  increase  in  transmitted  symbol  rate, 
producing  a  4  dB  decrease  in  processing  gain.  In  general  the  loss  in  processing 
gain  due  to  inclusion  of  a  guard  band  to  protect  against  an  offset  of  ±X%  of  the 
CVI  is 
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which  is  illustrated  in  Figure  9.1-1.  Of  course,  the  CVI  should  be  made  as 
small  as  possible  in  order  to  minimize  the  amount  of  time  the  jammer  has  to 
observe  and  process  a  particular  set  of  spreading  waveforms  and  then  use  the 
extracted  parameters  in  a  jamming  signal.  But  the  minimum  CVI  is  constrained 
by  the  amount  of  processing  gain  loss  which  can  be  tolerated  which  is  a  func¬ 
tion  of  both  the  CVI  length  and  the  maximum  timing  error  anticipated  in  the 
system.  A  major  disadvantage  of  the  hard  edge  CVI  approach  is  that  in  systems 
which  include  a  mix  of  high  precision  clock  units  and  low  precision  clock  units 
the  guard  band  is  usually  fixed  for  all  data  links  at  a  value  which  will  accom¬ 
modate  the  lowest  precision  clock.  Therefore  even  the  high  precision  clock 
units  pay  the  penalty  of  reduced  processing  gain.  This  situation  can  be 
avoided  through  a  modified  CVI  approach  in  which  changes  in  the  spreading 
waveform  set  are  evolutionary  rather  than  abrupt.  This  technique  is  referred 
to  as  soft  edge  CVI. 

9.2  Soft  Edge  CVI 

In  the  soft  edge  CVI  method  small  changes  are  continually  intro¬ 
duced  in  the  spreading  waveforms  throughout  the  CVI  in  such  a  way  that  the 
waveforms  completely  evolve  from  one  set  to  the  next  i">  CVI.  For  example  if 

in  a  hard  edge  CVI  system  the  "0"  symbol  is  represented  by  the  PN  segment  S  .i 
+■  h  Q 

during  the  i  CVI  interval  and  by  S  during  the  next  interval,  the  change 

0,1+1  th 

from  S  .  to  S  .+1  is  made  abruptly  at  the  end  of  the  i  interval.  In  an 

U  »  1  0  j  1  X 

equivalent  soft  edge  CVI  system  changes  are  made  one  chip  at  a  time  to  SQ  ■ 
so  that  at  the  end  of  the  CVI  s_  :  has  been  completely  changed  to  S  ...  With 

O  y 1  0)1  'X 

SQ  ^  and  SQ  orthogonal,  the  two  PN  segments  are  alike  in  50%  of  the  chip 
positions  and  different  in  50%  of  the  chip  positions.  Therefore  only  50%  of 
the  chips  need  be  changed  during  each  CVI.  and  they  can  be  selected  for  change 
in  a  pseudo-random  order  which  is  common  to  transmitter  and  receiver.  This 
evolutionary  change  therefore  provides  some  additional  processing  gain  against 
repeat  jammers  but  more  importantly,  it  provides,  (1)  a  graceful  degradation 
in  performance  as  a  function  of  timing  error  and,  (2)  performance  of  each  unit 
is  determined  by  the  clock  precision  of  the  unit,  not  the  precision  of  the 
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lowest  precision  clock  in  the  system.  These  advantages  are  illustrated  as 
follows.  Since  50%  of  the  chips  are  changed  in  a  CVI,  25':.  of  the  chips  have 
changed  by  the  mid-point  of  the  CVI.  Thus  75'  of  the  chips  in  the  mid-point 
"0"  symbol  match  the  initial  SQ  ■  PN  segment  and  likewise  75%  of  the  chips 
in  the  mid-point  "0"  symbol  match  the  final  SQ  .+1  PN  segment.  If  the  mid-point 
"0"  symbol  is  correlated  with  either  S  •  or  S  .,,,  75%  of  the  chips  agree, 

25%  disagree  so  that  the  correlation  value  is  one-half  the  maximum  vulue, 
which  indicates  a  loss  in  processing  gain  of  6  dB  for  a  unit  whose  timing  error 
is  50%  of  a  CVI.  For  a  unit  whose  timing  error  is  zero  the  mid-point  "0"  received 
symbol  would  be  correlated  with  the  same  PN  segment  so  that  there  is  no  loss  in 
processing  gain.  Thus  the  processing  gain  available  to  each  unit  is  a  function 
only  of  that  unit's  timing  accuracy.  Processing  gain  loss  as  a  function  of 
timing  error  is 

P91 oss *  dB  =  '20  log10  ■  T5o) 

where  X  is  timing  error  expressed  in  percent  of  CVI.  This  function  is  plotted 
in  Figure  9.2-1. 

Note  that  a  timing  offset  of  ±30%  in  the  soft  edge  CVI  system  results 
in  a  3.1  dB  loss  in  processing  gain  whereas  as  noted  earlier  a  guard  band  to 
protect  against  a  ±30%  offset  in  the  hard  edge  CVI  system  would  result  in  a 
4  dB  loss  in  processing  gain. 

9.3  Acquisition  Characteristics 

No  acquisition  or  reacquisition  operations  are  required  in  CVI  systems, 
and  hence  no  time  need  be  dedicated  to  acquisition,  whenever  the  relative  timing 
between  the  transmitter  and  receiver  is  within  one  CVI  for  hard  edge  systems  or 
within  some  fraction  of  a  CVI  for  soft  edge  systems.  Under  these  conditions 
the  MFs  in  the  receiver  are  matched  to  the  incoming  symbols  (or  nearly  matched 
in  the  case  of  soft  edge  systems)  and  symbol  detection  is  accomplished  immediately. 
Under  conditions  of  larger  timing  uncertainties  a  search  over  the  CVIs  in  the 
uncertainty  region  is  required.  But  note  that  in  DS  systems  the  acquisition 
search  is  conducted  in  4-chip  increments  and  the  timing  resolution  required 
for  data  demodulation  is  a  small  fraction  of  a  chip,  whereas  in  CVI  systems 
the  search  is  conducted  in  increments  of  some  fraction  of  a  CVI  and  timing 
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Figure  9.2-1.  Processing  Gain  Loss  in  Soft  Edge 

CVI  System 


need  be  resolved  to  only  some  fraction  of  a  CVI.  Since  a  CV1  is  typically 
many  orders  of  magnitude  greater  than  a  chip  width  the  search  time  is  many 
orders  of  magnitude  less  for  the  CVI  system.  And,  of  course,  the  need  for 
an  acquisition  search  is  greatly  reduced  in  CVI  systems  since  a  much  larger 
timing  error  can  be  tolerated. 
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This  investigation  has  focused  on  the  fundamental  principles  of 
the  acquisition  process,  mathematically  quantifying  potential  acquisition  per¬ 
formance  as  a  function  of  (1)  receiver  structures  and  devices,  (2)  spreading 
waveforms  and  (3)  specially  structured  acquisition  signals.  Certainly  the 
most  obvious  conclusion  is  that  MF  devices  offer  a  relatively  straightforward 
means  (for  DS  systems  at  least)  of  achieving  two  to  three  orders-of-magnitude 
improvement  in  acquisition  time.  These  magnitudes  of  improvement  are  quite 
significant  for  reacquisition  operations  in  which  only  "coarse"  or  "fine" 
timing  uncertainties  exist.  In  addition,  acquisition  with  a  MF  results  in 
the  smallest  residual  timing  uncertainty  at  the  completion  of  the  acquisition 
process.  Therefore,  every  attempt  should  be  made  to  exploit  this  rapidly 
developing  technology,  especially  for  application  to  DS  data  links  which 
require  both  rapid  reacquisition  and  precision  ranging  capabilities.  A  spec¬ 
ific  recommendation  for  further  work  in  this  area  is  presented  in  the  next 
section. 

Insofar  as  the  influence  of  the  spreading  waveform  is  concerned 
direct  search  of  FH  was  shown  to  provide  the  same  acquisition  speed  advan-  - 

tage  over  direct  search  of  DS  as  did  a  MF  acquisition  of  DS.  However,  it  was 

also  seen  that  direct  search  of  FH  results  in  a  larger  residual  timing  uncer¬ 
tainty  at  the  completion  of  the  acquisition  process  than  does  a  MF  acquisition 
of  DS.  Basically  direct  search  of  FH  has  an  acquisition  speed  advantage  over 
direct  search  of  DS,  for  the  same  initial  timing  uncertainties,  simply  because 
the  timing  accuracy  required  for  data  demodulation  in  FH  systems  is  less  than 

that  required  in  DS  systems.  If  the  only  purpose  of  the  spread  spectrum  link 

is  to  provide  data  transfer,  then  a  FH  system  employing  direct  search  acquisition 
has  an  acquisition  speed  advantage  over  a  DS  system  employing  direct  search 
acquisition.  If  precision  ranging  is  a  requirement  then  the  FH  system  must 
employ  a  MF  structure  matched  to  the  frequency  hopping  pattern  which  will  re¬ 
sult  in  acquisition  time  and  residual  timing  uncertainty  comparable  to  that 
provided  by  a  DS  system  employing  a  MF  receiver.  The  choice  between  FH  and 
DS  then  must  take  into  consideration  the  relative  complexity  of  the  two  MF 
receiver  structures,  and  here  the  DS  system  has  a  decided  advantage.  Evaluation 
of  the  FH/DS  waveform  for  acquisition  indicated  that  its  performance  falls 
between  that  of  DS  and  FH  with  regard  to  both  acquisition  time  and  residual 
timing  uncertainty.  In  addition,  its  receiver  structure  is  more  complex  than 
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either  of  the  others  so  that  it  cannot  be  considered  an  attractive  approach 
to  acquisition.  The  bottom  Tine  on  these  comparisons  is  that  essentially  the 
same  acquisition  performance  can  be  achieved  from  DS  and  FH  waveforms  if  the 
requirement  is  to  resolve  timing  from  the  same  initial  value  to  the  same  final 
value;  and  FH/DS  can  be  used  to  achieve  similar  acquisition  speed  but  with 
larger  residual  timing  uncertainty.  Consequently,  the  choice  of  spreading 
function  for  a  particular  data  link  application  should  be  much  more  strongly 
influenced  by  other  system  requirements  than  by  acquisition  performance. 

Several  approaches  to  specially  structured  acquisition  signals 
are  presented  in  the  report.  These  methods  are  of  particular  interest  for 
initial  acquisition  or  for  reacquisition  with  relatively  large  timing  uncer¬ 
tainty,  that  is,  conditions  under  which  even  the  acquisition  performance  of  a 
MF  is  not  sufficient.  As  often  noted  in  discussions  of  these  methods,  any 
structure  in  the  transmitted  signal  can  be  used  by  a  sufficiently  intelligent 
jammer  to  increase  his  effectiveness.  These  techniques  should,  therefore, 
only  be  considered  if  acquisition  or  reacquisition  with  total  or  very  large 
timing  uncertainty  is  a  real  system  requirement,  and  if  the  potential  vul¬ 
nerability  is  acceptable  in  light  of  the  anticipated  ECM  environment.  Given 
that  the  requirement  exists  and  that  the  potential  vulnerability  is  acceptable, 
the  "Combined  Search  and  Data  Transmission  Method"  presented  in  Sections  8.3 
through  8.5  provides  a  straightforward  means  of  trading  between  waveform 
structure  and  acquisition  time.  In  this  method  a  short  code  spread  signal 
if  first  acquired  (preferably  by  a  MF)  and  then  the  long  code  seed  is  trans¬ 
mitted  as  data  on  the  short  code  spread  signal.  As  shown  in  Section  8.5,  the 
length  of  the  short  code  can  be  selected  to  minimize  the  acquisition  time,  with 
longer  short  code  lengths  resulting  in  greatly  increased  acquisition  time. 

An  alternative  to  conventional  short  codes  are  the  component  codes  described 
in  Appendix  B.  These  codes  are  generated  by  combining  several  codes  in  such 
a  way  that  the  period  of  the  resulting  code  is  equal  to  the  least  common 
multiple  of  the  periods  of  all  the  components  but  that  the  acquisition  time 
is  essentially  the  time  to  acquire  the  longest  component  only.  The  result  is 
an  acquisition  time  significantly  less  than  that  of  a  comparable  length  PN 
code  but  the  acquisition  detector  requires  an  acquisition  device  for  each  code 
component  in  order  to  fully  realize  the  potential  speed  advantage.  Further 
investigation  of  the  implementation  and  vulnerability  of  these  codes  is  required 
before  a  meaningful  comparison  can  be  made  between  them  and  conventional  short 
code  techniques. 
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The  special  waveform  techniques  which  have  the  potential  for  making 
the  greatest  impact  on  tactical  AJ  data  links  are  the  CV I  techniques  discussed 
in  Section  9.  These  techniques  allow  data  demodulation  with  relatively  large 
timing  uncertainty  between  transmitter  and  receiver,  thereby  avoiding  any 
reacquisition  operations  for  most  tactical  scenarios.  Direct  utilization  of 
the  CVI  methods  would  make  the  data  link  somewhat  more  susceptible  to  smart 
jamming  but  the  inclusion  of  well  designed  message  validation  techniques  can 
greatly  reduce  this  vulnerability.  These  are  areas  for  further  study  and  are 
discussed  in  the  next  section. 
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11.0  RECOMMENDATIONS  FOR  FUTURE  WORK 

Having  well  established  the  fundamental  principles  of  acquisition 
and  identified  areas  of  potential  performance  advantages,  it  is  appropriate  to 
now  exploit  the  most  promising  techniques.  The  "ultimate"  acquisition  strategy 
is  one  which  is  fully  AJ  protected  at  all  times  and  which  provides  instantaneous 
acquisition.  Such  a  strategy  would  be  the  basis  for  a  highly  robust  model  which 
could  meet  a  wide  range  of  tactical  communication  requirements,  including  "push- 
to-talk"  requirements.  Such  an  acquisition  strategy  is  obviously  sometime  off 
in  the  fucure  but  does  represent  a  most  significant  and  an  achievable  goal,  and 
this  section  outlines  a  plan  for  achieving  that  goal. 

Matched  filter  technology  is  recognized  as  one  key  to  improved  modem 
performance,  for  both  acquisition  and  data  demodulation.  Devices  which  are 
either  presently  available  or  under  development,  CCDs,  digital  correlators,  SAW 
convolvers  and  the  hybrid  array  correlator,  were  described  in  this  report  and 
their  utilization  and  performance  as  acquisition  processors  indicated.  Inves¬ 
tigation  of  new  and  more  powerful  MF  concepts,  such  as  the  optical  processing 
techniques  described  in  Appendix  D,  should  be  encouraged  and  supported  by  the 
AJ  systems  community.  The  benefit  to  the  AO  systems  community  would  not  only 
be  the  continual  evaluation  of  the  applicability  of  these  new  devices  to  AJ 
systems  but  also  that  development  of  the  most  promising  concepts  would  be 
driven  such  that  the  end  product  is  directly  applicable  to  future  AJ  systems. 

These  devices  would  represent  powerful  signal  processing  resources  which  could 
support  multiple  receiver  functions  including  acquisition,  data  demodulation, 
carrier  tracking  and  estimation  of  other  pertinent  system  parameters. 

Concurrent  with  the  investigation  of  new  MF  devices,  the  performance 
advantages  of  present  state-of-the-art  devices  should  be  verified  and  demonstrated 
through  implementation,  at  the  earliest  opportunity,  in  a  state-of-the-art  AJ 
system.  The  Army's  Modular  Integrated  Communication  and  Navigation  System  (MICNS) 
presently  in  development  at  Harris  represents  an  excellent  opportunity  for  near 
term  application  of  a  state-of-the-art  MF  to  acquisition  in  a  state-of-the-art 
DS  system.  Demonstration  with  the  MICNS  would  be  a  relatively  quick  turnaround 
effort  requiring  a  minimal  design  and  breadboard  effort  for  peripheral  subsystems. 
This  application  would  then  provide  the  design  experience  and  test  data  upon  which 
to  establish  the  requirements  for  new  MF  devices.  This  implementation  activity 
would  be  supplemented  by  an  analysis  effort  aimed  at  optimizing  the  effective¬ 
ness  of  the  MF  based  acquisition  processor.  The  correlation  period  of  MFs  is 
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limited  by  physical  constraints  on  the  devices  and  by  signal  decorrelation 
due  to  Doppler.  In  most  applications  the  correlation  period  will  be  between 
10  and  20  microseconds.  For  low  symbol  rates  more  than  one_xorrelation  per 
acquisition  test  will  be  required  in  order  to  obtain  a  sufficiently  high  prob¬ 
ability  of  acquisition  detection  at  typical  jammer  to  signal  ratios.  Therefore, 
the  acquisition  processor  must  include  circuitry  for  sampling  and  processing 
the  MF  output.  The  process  of  converting  the  MF  output  to  digital  samples 
involves  trade-offs  between  sample  rate  and  sample  quality,  pulse  shaping 
techniques,  AGC,  etc.  The  processing  of  the  digital  samples  involves  trade¬ 
offs  between  processing  gain,  probability  of  detection  and  probability  of 
false  alarm,  and  complexity  of  the  processing  hardware.  Recent  advances  in 
LSI  technology  permit  quite  sophisticated  processing.  For  example,  it  may 
be  practical  to  keep  track  of  past  false  alarm  patterns  and  to  use  this  infor¬ 
mation  to  dramatically  reduce  the  dismissal  time  of  future  false  alarms.  In 
addition,  the  processor  could  be  adaptive  or  programable  such  that  its  para¬ 
meters  (processing  gain,  sampling  rate,  etc.)  could  be  changed  to  be  compatible 
with  the  environment,  either  automatically  or  by  user  control.  The  inclusion 
of  these  signal  processing  capabilities  into  a  MF  acquisition  processor  would 
result  in  a  greatly  enhanced  acquisition  strategy,  and  could  be  accomplished 
with  today's  technology. 

The  next  step  towards  the  "ultimate"  acquisition  strategy  is 
the  formulation  of  receiver  structures  and  waveform  designs  which  will  exploit 
the  high  speed,  high  density  digital  processing  capabilities  of  the  forthcoming 
VHSIC  and  VLSI  technologies.  These  techniques  will  provide  maximum  processing 
power  at  minimum  cost  by  performing  most  of  the  receiver  functions  in  the 
relatively  low  cost  digital  domain  rather  than  in  the  relatively  high  cost 
analog  domain.  The  high  speed  of  the  devices  will  permit  considerable  pro¬ 
cessing  of  the  received  signal  for  enhanced  message  identification  and  smart 
jammer  dismissal.  This  processing  power  will  allow  the  use  of  more  complex 
acquisition  codes  which  will  give  a  decided  edge  to  the  friendly  receiver  over 
the  intelligent  jammer.  "Long"  short  codes  and  more  complex  component  codes 
will  be  more  attractive  for  acquisition.  But  the  most  attractive  signalling 
techniques  for  this  era  of  technology  appear  to  be  the  CVI  techniques.  The 
basic  techniques  rely  on  a  sufficiently  high  rate  of  code  change  to  reduce 
susceptibility  to  smart  jammers.  Shortening  the  CVI  increases  the  protection 
against  smart  jammers.  However,  clock  instability  and  vehicle  dynamics  place 
a  limit  on  how  short  the  CVI  can  be  made  and  still  maintain  the  desired  processing 


i 


11-2 


gain.  If  moderately  accurate  position  information  is  available,  the  clock  is 
the  limiting  factor  and  the  length  of  the  CVI  is  a  trade-off  between  clock 
accuracy  and  speed  of  the  smart  jammer.  However,  with  the  processing  power 
of  VHSIC  and  VLSI  message  validation  techniques  can  be  incorporated  to  signi¬ 
ficantly  reduce  the  vulnerability  of  a  relatively  long  CVI. 

The  final  step  in  the  plan  is  the  merging  of  high  accuracy  clocks, 
high  precision  position  estimation  and  powerful  digital  signal  processing. 

At  that  point  the  CVI  can  be  reduced  and  the  small  CVI,  supported  by  the  power¬ 
ful  message  validation  procedures ,  will  result  in  the  "ultimate  acquisition 
strategy  -—  full  AJ  protection  and  instantaneous  acquisition. 

Harris  is  prepared  to  discuss,  in  detail,  >tudy  activities  which 
will  support  the  areas  outlined  above. 
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FH  ACQUISITION  BY  ASEAT 


An  interesting  technique  for  acquisition  of  FH,  the  Autoregressi ve 
Spectral  Estimation  Acquisition  Technique  (ASEAT),  has  recently  been  reported 
in  the  literature  [l,  2].  A  simplified  block  diagram  of  the  system  is  shown 
in  the  attached  figure  with  operation  as  follows.  The  N-stage  linear  feed¬ 
back  shift  register  (LFSR)  in  the  transmitter  provides  the  N-bit  address  for 

the  frequency  to  be  transmitted  during  each  frequency  interval.  The  allows 
N 

2  different  frequencies  for  hopping.  Let  the  contents  of  the  LFSR  during 
the  ktb  interval  be  denoted  by 


bN  bN-l 


b 


1 


where  hj  ,s  the  most  significant  bit  (MSB)  and  b„  is  the  least  significant 
(LSB).  The  corresponding  transmitted  frequency  during  this  interval  is  f 
the  (HI)  interval  the  LFSR  is  shifted  one  bit  so  that  the  contents  are* 
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bN+l  bN  b2 

producing  the  frequency  fk+J.  For  transmission  intervals  k  through  k+N-1  the 
LFSR  contents  and  corresponding  frequencies  are: 
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Acquisition  by  ASEAT 


Note  that  the  list  of  MSBs  constitute  the  LFSR  contents  during  the  initial 
interval.  Thus,  if  the  receiver  can  estimate  the  received  frequencies  f^, 

fk+l»  - -  ^k+N-i  with  sufficient  accuracy  to  have  the  MSB  of  the  addresses 

correct,  then  the  receiver  will  exactly  know  the  initial  contents  of  the 
transmitter' s  LFSR.  Since  all  future  contents  of  a  LFSR  are  uniquely  deter¬ 
mined  from  the  contents  at  any  one  instance,  the  receiver's  LFSR  can  be 
loaded  so  as  to  match  that  of  the  transmitter.  Thus  acquisition  is  achieved 
in  N  intervals,  by  sending  one  bit  of  timing  information  (b^  then  b2  then  b^ 

....  then  bN)  during  each  interval.  If  the  LFSR  is  connected  to  produce  a 
w  N 

maximal -length  sequence  the  length  is  2  -  I.  Therefore,  this  technique  is 

Ml 

seen  to  be  a  data  transmission  method  in  which  acquisition  of  a  2  chip 
long  code  is  achieved  by  transmission  of  N  bits  of  timing  data.  Unfortunately 
a  smart  jammer  can  also  use  the  ASEAT  to  obtain  synchronization  with  the  trans¬ 
mitter  and  can  then  exactly  produce  the  FH  pattern  with  any  desired  delay  or 
advance.  Thus  such  a  system  has  a  rapid  acquisition  capability  but  is  also 
quite  vulnerable  to  smart  jamming.  Normally  the  technique  used  to  overcome 
this  torn  of  vulnerability  is  nonlinear  code  generation,  in  which  case  the 
MSBs  of  the  sequence  of  received  frequencies  will  not  be  the  initial  contents 
of  the  transmitter's  shift  register.  Consequently  the  ASEAT  does  not  appear 
to  be  a  viable  acquisition  method  for  AJ  systems. 
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Component  codes  offer  a  means  of  introducing  structure  into  a 
DS  spreading  signal  in  such  a  way  that  it  can  be  exploited  by  the  receiver 
to  reduce  acquisition  time.  As  emphasized  throughout  the  report,  any  reduc¬ 
tion  in  randomness  of  the  spreading  waveform  increases  the  vulnerability  to 
smart  jammers  so  that  the  use  of  component  codes  does  reduce  the  AJ  protection. 
On  the  otV'r  hand  they  represent  an  alternative  to  the  other  structured  wave¬ 
form  methods,  such  as  conventional  short  codes,  and  may  be  employed  when 
acquisition  is  to  be  performed  in  an  environment  of  relatively  unsophisticated 
jammers . 

A  component  code  (S)  is  one  that  is  formed  by  performing  some 
Boolean  function  (B)  on  several  component  codes  (X.- ),  that  is 

S  =  B  (Xj,  X2,  ....  Xn) 

The  period  of  the  code  S  is  equal  to  the  least  common  multiple  of  the  periods 
of  all  component  codes.  If  the  component  code  periods  are  chosen  to  be 
relatively  prime,  the  period  of  S  is  equal  to  the  product  of  the  periods  of 
the  component  codes.  Furthermore,  if  B  is  chosen  correctly  the  code  S  will 
have  nearly  equal  numbers  of  1 1 s  and  O's  and  the  component  code  will  exhibit 
a  detectable  cross  correlation  with  each  of  its  components.  It  is  this  last 
property  which  allows  the  component  code  to  speed  up  the  acquisition  process. 

At  this  point  an  example  is  in  order.  Consider  a  three  component 
code  and  the  Boolean  function 

S  =  Xj  x  X2  +  X2  x  X3  +  X3  x  Xj 

where  Xj,  X2  and  Xg  are  the  component  codes.  We  could  choose  the  component 
code  periods  as 

Pv  =  27  -  1  =  127  (prime) 

X1 

Pv  =  29  -  1  =  511  (*  7  x  73) 

X2 

Pv  =  210  -  1  =  1023  (=  3  x  11  x  31). 

*3 


Since  these  periods  have  no  common  factors  they  are  relatively  prime,  thus  the 
length  of  the  composite  code  is 


P$  =  127  x  511  x  1023  =  66389631. 

This  period  is  somewhat  shorter  than  the  period  of  a  maximal  length  sequence 
having  7  +  9  +  10  =  26  state  variables,  but  longer  than  that  of  a  25  state 
variable  maximmal  length  sequence.  For  this  particular  code  there  are  358785 
more  l's  than  0's.  It  would  appear  that  a  maximal  length  sequence  with  26 
state  variables  is  preferable  to  this  composite  code  on  the  basis  of  slightly 
longer  period  and  better  1/0  balance.  This  is  true  until  one  considers  the 
acquisition  process.  Using  the  proper  Boolean  function  allows  acquisition  of 
each  of  the  component  codes  independently  and  simultaneously.  Consider  what 
happens  when  the  composite  code  (transmitted  as  1— ►+  1  and  0— ^-l)  is 
correlated  with  a  replica  of  any  of  its  components.  This  gives 

R£s  =  2  E[X©S]  -  1 

due  to  the  translation  0 — ►  -!.  If  X  is  a  time  shifted  version  of  a  component 
(essentially  uncorrelated)  then 

RXS  ~  0 

while  if  X  is  identical  to  a  component 

Rxs  ~ 

Thus  the  problem  of  acquiring  one  component  is  simply  the  problem  of  acquiring 

a  short  code  only  with  6  dB  less  signal  power  (neglecting  code  self  noise). 

Since  all  components  may  be  acquired  in  parallel,  the  acquisition  time  is 

essentially  equal  to  the  time  required  to  acquire  the  longest  components.  For 

our  example  code,  the  number  of  (half  chip)  cells  to  search  is  reduced  from 
27  11 

~2  down  to  ~2  for  a  reduction  of 
~216  =  65536. 

Of  course  acquisition  is  not  65000  times  faster  since  the  cells  are  searched 
with  a  6  dB  lower  S/0,  but  considerably  accelerated  acquisition  is  possible. 
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TIMING  UNCERTAINTY  AT  REACQUISITION 


The  time  required  to  reacquire  a  spread  spectrum  link  after  a  period 
of  link  outage  is  directly  proportional  to  the  timing  uncertainty  which  has 
accrued  during  the  outage.  It  is,  therefore,  of  interest  to  identify  the 
sources  of  timing  error  and  their  relative  contributions  to  the  total  error. 

In  addition,  an  estimate  of  the  total  timing  error  is  required  to  properly 
initialize  the  reacquisition  processor.  In  a  direct  search  reacquisition 
system  the  accrued  timing  error  defines  the  uncertainty  region  which  must  be 
searched;  in  a  MF  reacquisition  system  it  dictates  the  amount  by  which  the 
local  reference  must  be  advanced.  Underestimating  the  timing  error  may  result 
in  a  missed  acquisition;  overestimating  may  result  in  unnecessarily  long  reacq¬ 
uisition  time.  The  total  timing  uncertainty,  T  ,  is  a  function  of  PN  clock 
offset,  vehicle  dynamics,  error  growth  rate  of  any  inertial  navigation  system 
(INS)  used  for  position  estimation  and  loss-of-link  time. 

The  total  timing  error  that  accrues  between  two  previously  synchron¬ 
ized  clocks  is  due  to  clock  offset,  aging  and  a  randomly  varying  term  which  is 
characterized  by  the  Allan  variance.  Clock  offset  consists  of  initial  setting 
error,  temperature  stability,  etc.  and  produces  a  timing  error  which  is  a  linear 
function  of  loss-of-link  time.  Error  due  to  aging  can  usually  be  neglected  for 
the  relatively  short  periods  of  link  outage  encountered  in  tactical  communi¬ 
cation  systems.  The  randomly  varying  term  is  also  usually  negligible  in  com¬ 
parison  to  the  offset  error  for  clocks  of  the  quality  used  in  tactical  systems. 
Thus,  the  timing  uncertainty  contributed  by  the  clocks  is  essentially  that 
produced  by  clock  offset  and  is  a  linear  function  of  loss-of-link  time.  For 
example,  considering  the  transmitter's  clock  to  be  a  perfect  reference,  a 
clock  offset  in  the  receiver  of  10"^  would  produce  1  microsecond  of  timing 
offset  per  second  of  link  outage.  For  MF  reacquisition  this  would  indicate 
that  the  local  reference  should  be  advanced  1  microsecond  for  each  second  of 
link  outage,  whereas  for  a  direct  search  system  the  uncertainty  region  would 
grow  at  a  rate  of  2  microseconds  per  second  since  the  direction  of  relative 
offset  is  unknown.  Note  that  for  this  source  of  timing  error,  as  is  true  for 
the  other  sources  as  well,  the  uncertainty  region  for  direct  search  is  twice 
as  large  as  the  reference  advance  for  the  MF.  However,  the  MF  may  have  to 
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wait  an  amount  of  time  equal  to  twice  its  advance  before  coincidence  occurs. 
In  the  following  discussions  of  timing  uncertainties  only  the  one-sided  error 
will  be  indicated,  i.e.,  the  amount  of  MF  advance. 

As  noted  above  the  timing  uncertainty  due  to  the  PN  clock  is  a  linear 
function  of  the  loss-of-link  time 


uc 


"  Co  T1  x  10 


+6 
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where  CQ  is  the  clock  offset,  T-j  is  the  loss-of-link  time  in  seconds  and  Tuc 
is  the  timing  uncertainty  produced  by  the  clock  in  microseconds.  The  timing 
uncertainty  due  to  differential  velocity  between  transmitter  and  receiver  is 
also  a  linear  function  of  loss-of-link  time,  with  a  rate  of  approximately  1.1 
microseconds  per  second  at  Mach  1.  Therefore,  with  differential  velocity  v 
in  Km/hr  the  timing  uncertainty  in  microseconds  is 


uv 


=  3.34  x  10‘3  v  T 


1 
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However,  a  second  order  tracking  loop  can  accurately  track  out  constant  velo¬ 
city,  even  over  relatively  long  periods  of  link  outage,  so  that  this  error 
may  not  be  a  significant  factor.  Differential  acceleration  produces  timing 
uncertainty  which  is  porportional  to  the  square  of  T^.  For  a  differential 


acceleration  of  "a"  gs. 


ua 


=  32.0  x  10“3  a  T}2 
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where  again  the  timing  uncertainty  is  in  microseconds. 

The  timing  uncertainties  due  to  clock,  differential  velocity  and 
differential  acceleration  are  compared  in  Figure  B-l  where  the  assumed  system 
parameters  are: 

1-7 


=  10 


200  Km/hr 
2  gs 


It  is  noted  that  velocity  is  the  driver  for  loss-of-link  times  up  to  about 
11  seconds.  However,  as  noted  earlier,  constant  velocity  can  be  tracked  out 
by  a  second  order  loop  while  the  link  is  established  and  very  accurately  esti¬ 
mated  while  the  link  is  lost.  Thus,  the  major  contributions  are  actually 
acceleration  and  clock  offset.  Of  these  two,  clock  offset  is  the  major  factor 
for  the  first  few  seconds  and  then  acceleration  becomes  the  major  factor. 
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Timing  uncertainty  due  to  acceleration  can  be  greatly  reduced  by  INS 
aided  position  estimation.  A  relatively  low  cost  INS  unit  typically  employed 
in  RPVs  can  provide  range  estimation  with  a  one  a  error  of 
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6  x  10'3  T-j3 
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where  are  is  in  feet  and  T-j  is  in  seconds.  This  produces  a  PN  code  timing 
uncertainty  (+3a)  of 


uINS 


=  18  x  10"6  T^ 
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where  T^j^  is  in  microseconds  and  T^  in  seconds.  With  such  a  unit  providing 
range  information  the  timing  uncertainties  due  to  velocity  and  acceleration 
are  replaced  by  T  and  the  total  timing  uncertainty  is  due  to  the  clock 
and  the  INS.  Obviously  the  clock  offset  is  the  dominant  factor  for  most  link 
outages.  Equating  Tuc  and  TuI^  and  solving  for  T-j  indicates  that  the  INS 
does  not  become  the  dominant  factor  until  the  link  has  been  out  for  over  1.5 
hours.  Thus,  in  the  case  of  INS  aided  positioning  the  clock  offset  is,  in 
reality,  the  only  significant  factor  in  timing  uncertainty  for  reasonably  long 
link  outages. 

The  maximum  acceptable  timing  uncertainty  at  reacquisition  is  an 
important  system  parameter  in  that  it  establishes  the  minimum  update  rate. 

If,  for  example,  the  maximum  uncertainty  is  10  microseconds  and  the  total 
uncertainty  is  due  to  clock  offset  and  differential  acceleration  then  the  max¬ 
imum  period  of  link  outage  is  11.7  seconds  for  the  example  cited  above.  If 
INS  aided  positioning  is  available  only  clock  offset  need  be  considered  and 
the  maximum  T-j  becomes  100  seconds.  On  the  other  hand,  if  a  CVI  technique  is 
employed  with  a  CVI  length  of  hundreds  of  milliseconds  the  system  can  operate 
for  days  between  timing  updates  based  on  these  considerations. 
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In  addition  to  the  CCD  and  SAW  convolvers,  digital  correlators 
and  hybrid  array  technologies  cited  in  Section  5,  optical  approaches  warrant 
further  study  for  performing  high-speed,  parallel  processing.  The  key  element 
in  optical  processing  is  a  Bragg  cell  which,  in  most  implementations  functions 
as  an  analog  or  a  digital  delay  line.  Generally  these  devices  have  bandwidths 
in  the  50-300  MHz  range,  delay  times  ranging  from  1  to  40  sec  and  time-band¬ 
width  products  in  the  range  of  1000-2000. 

For  some  communication  applications  the  integration  time,  which  is 
basically  the  delay  time,  is  too  short.  In  recent  years,  however,  a  new  class 
of  processing  architectures  have  evolved  in  which  the  integration  time  can  be 
extended  appreciably  to  provide  additional  correlation  gain.  This  class  of 
optical  systems  are  generally  called  time-integrating  optical  systems  (TIC). 

One  attractive  feature  of  optical  processing  is  that  the  operations  are  naturally 
performed  in  parallel.  Of  equal  significance  for  rapid  acquisition  is  the  fact 
that  optical  systems  are  two-dimensional  so  that  the  orthogonal  axis  can  be  used 
to  display  doppler  if  the  principal  axis  is  used  to  display  time  delay. 

We  first  describe  an  optical  system  that  performs  a  matched  filtering 
operation  similar  to  that  described  before;  this  example  will  illustrate  some  of 
the  basic  principles  of  optical  processing.  Figure  D-l  shows  an  optical  system 
in  which  light  from  a  point  source  is  collimated  by  lens  L^.  The  source  might  be 
a  small  solid  state  light  emitting  diode  or  an  injection  laser  diode.  Suppose 
that  we  modulate  the  source  with  the  PN  sequence  f1(t)=C(t)  and  that  this  light 
illuminates  a  Bragg  cell  which  is  driven  by  the  received  signal  f2(t)  which  is 
the  same  sequence,  with  a  large  starting  time  uncertainty,  plus  noise.  For  the 
moment,  suppose  we  treat  the  noise  free  case;  the  Bragg  cell  operates  so  that 
the  waveform  within  the  cell  is  C(t-x/v)  where  x  is  a  space  coordinate  in  the 
plane  of  the  Bragg  cell  and  v  is  the  velocity  of  the  acoustic  wave.  The  light 
emerging  from  the  Bragg  cell  is  then  C(t)C(t+x)  where  x=-x/v. 

If  we  image  the  Bragg  cell  onto  a  linear  photodetectora  array  contain¬ 
ing  N  elements,  we  can  associate  a  given  time  delay  .  with  the  jth  element  in 
the  array.  If  we  integrate  for  a  time  T,  where  T  may  be  much  longer  than  the 
delay  time  of  the  Bragg  cell,  we  have  that  the  output  of  the  jth  element  is 

T 

Rj  «  R(tj)  =  f  C ( t ) C ( t+r j )dt  . 

*'o 
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Since  the  number  of  elements  in  the  array  may  be  of  the  order  of  1000  and  if 
we  measure  the  correlation  at  h  chip  intervals,  we  see  that  we  can  dismiss  500 
chip  intervals  in  each  time  interval  of  duration  T. 

By  using  multiple  sources  or  multi-channel  Bragg  cells,  along  with  a 
multi-channel  or  two  dimensional  photodetector  array,  we  can  increase  the  time 
delay  searched  in  parallel  by  a  corresponding  factor.  Alternatively,  different 
dopplers  could  be  provided  in  each  channel  with  the  same  PN  sequence  phasing  to 
give  a  simultaneous  doppler/delay  output. 

It  should  be  noted  that  this  optical  system  performs  the  matched 
filtering  operation  without  the  need  to  time  reverse  the  reference  PN  sequence 
which  is  typical  of  many  other  approaches.  Very  long  sequences  can  therefore 
be  correlated  easily. 

If  it  is  desirable  to  rapidly  acquire  a  PN  sequence  when  the  doppler 
range  is  large  (or  when  fine  doppler  increments  must  be  tested),  a  different 
optical  configuration  might  be  useful.  Such  a  system,  not  shown  here,  is  called 
a  "triple-product  processor."  Its  configuration  is  similar  to  that  shown  in 
Figure  X  except  that  the  first  Bragg  cell  is  imaged  onto  a  second  Bragg  cell 
that  is  oriented  along  the  axis  orthogonal  to  the  first  cell.  Both  cells  are 
then  imaged  onto  a  two-dimensional  photodetector  array.  This  system  has  the 
interesting  property  that  it  displays  the  ambiguity  function  as  given  by  the 
combination  of  equations  (3.4-8)  and  (3.4-9).  The  integrated  output  is  given 
as 


R(T,(Od)  = 


I 

/ 


C(t)  C(t+r)e  Q  dt 


where  t  is  a  continuous  variable  in  the  horizontal  direction  and  is  a  con¬ 
tinuous  variable  in  the  vertical  direction.  For  a  1000x1000  photodetector  array 
we  can  examine  of  the  order  of  (500)  =250,000  doppler/delay  bins  in  a  time  inter 
val  of  T. 

Other  optical  configurations  can  be  used  to  handle  frequency  hopped 
signals.  For  example,  it  is  well-known  that  a  coherently  illuminated  optical 
system  displays  the  instantaneous  Fourier  transform  of  that  portion  of  a  signal 
resident  within  the  Bragg  cell  at  any  instance  in  time.  As  the  frequency  is 
hopped,  the  position  of  the  focused  beam  in  the  Fourier  plane  also  hops,  with  a 


D-2 


unique  one-to-one  mapping  of  hop  frequency  and  spatial  position.  Schemes  for 
activating  photodetectors  based  on  the  hopping  sequence  or  the  use  of  a  second 
Bragg  cell  to  introduce  a  reference  of  the  hopped  sequence  may  be  useful  to 
dehop  and  despread  the  signal  in  a  single  system. 

To  date  very  little  effort  has  been  applied  toward  investigating 
the  potential  uses  of  the  features  of  two-dimensional  or  multi-channel  optical 
processing.  We  propose  to  examine  this  potential  by  studying  and  adding  to  the 
existing  class  of  optical  processing  architectures  for  implementing  rapid  acqui¬ 
sition  algorithms.  The  focus  will  be  to  develop  techniques  that  provide  high¬ 
speed,  parallel  searches  for  doppler  and  time  delay  simultaneously.  Although 
we  currently  envision  using  the  optical  techniques  primarily  for  acquisition, 
there  are  also  possibilities  for  demodulation  and  detection  as  well.  Efforts 
will  also  be  directed  toward  system  configurations  that  are  compact  and  low-cost, 
but  these  considerations  are  of  less  importance  at  this  time. 


Figure  0-1.  Time-Integrating  Correlator 
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